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Thermal Degradation of Polymers At Low Rates 
Samuel L. Madorsky 
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In measuring the rates of thermal degradation of polymers at constant temperatures 
in the range of 200° to 500° C, one encounters a distortion of the initial parts of the rate 
curves during the time required to heat the specimen to the pyrolysis temperature. This 
distortion has been largely overcome by using a new apparatus, based on an electronic 
balance, in which the degradation of the specimen is carried out at very low rates. The 
apparatus is equipped with an electronic thermostat that controls the temperature to within 

0.1° C, and also with an automatie recorder that measures both the temperature and the 
loss of weight of the specimen. Pyrolysis experiments were carried out in a vacuum in 
this apparatus on pure specimens of polystyrene, poly(a-methylstyrene), polymethylene, 
high-molecular-weight polyethylene, and low-molecular-weight polyethylene. The initial 
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rates of degradation, in percent of sample per minute, varied from 0.008 to 0.095. 
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for the initial parts of the rate curves, in some eases up to about 30 percent loss of weight, 


the rate curves 
degradation. 


1. Introduction | 


| 

Studies on rates of thermal degradation of polv- | 
mers in a vacuum by the loss-of-weight method have | 
been carried out at the National Bureau of Standards 
for anumber of vears. The method [1, 2]? consisted 
in heating a 4 to 6 mg sample of polymer in a plati- 
num crucible suspended from a tungsten spring 
balance [1,5] and observing at intervals the loss of 
weight sustained by the sample, as indicated by the 
position of a crossline on the extension of the spring. 
For convenience, the rate experiments were usually 
planned to be completed within about 8 to 9 hrs. 
This placed a lower limit on the temperature of 
operation. In the range 250° to 500° © used in 
these studies it usually took about 15 min to heat 
the polymer sample from room temperature to that 
required in a given experiment. Some vaporization 
of the sample, in some cases [2] as much as 30 percent 
of its weight, took place during this heating-up 
period. Constancy of temperature in any given 
experiment was a required condition in all the rate 
studies. The initial loss thus caused a distortion of 
the initial part of the rate curve and masked the 
true picture of the initiation of thermal degradation. 

This difficulty could be largely overcome by 
carrving out the rate experiments for longer periods 
of time at lower temperatures so that the initial loss 
of sample would be reduced to as low a percentage 
as practicable. However, continuous operation for 
long periods of time suggests the use of automatic 
recording of loss of weight of sample. The spring 
balance did not lend itself easily automatic 
recording. 

A new apparatus for the study of rates of thermal 
degradation of polymers has been designed and con- 
structed. It records continuously and automati- 
cally time, temperature, and weight of the polymer 
sample undergoing decomposition. This new appa- 
ratus was used in a study of thermal degradation at 
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? Figures in brackets indicate the literature references at the end of this paper. 


at low rates resemble closely those previously obtained at high rates of 


low rates of several polymers. The polymers 
selected for this study were: Polystyrene, which 
vields both monomer and larger molecular frag- 
ments; polv(a-methylstvrene), which vields monomer 
exclusively; and nonbranched and branched poly- 
ethvlenes, which vield a spectrum of fragments 
containing from 1 to n carbons. 


2. Apparatus and Experimental Procedure 


The elements of the new apparatus (exclusive of 
the vacuum system) are shown diagramatically in 
figure 1, and the whole assembly is shown photo- 
graphically in figure 2.8. The equipment is_ built 
around an electronic microbalance, which is shown 
in figure 1 in two parts: A, the balance; and B, the 
electronic unit. From the arm P of the balance a 
crucible Q containing the polymer sample is sus- 
pended within the lower furnace F. The funetion 
of the upper furnace FE is deseribed later in this 
paper. The output of the recording balance is fed 
into one channel of a three-channel recording po- 
tentiometer C. The temperature of the furnace F 
is controlled to about +0.1° © by an electronic 
thermostat J, which senses the furnace temperature 
by means of a platinum resistance thermometer T 
built into the furnace. 

The temperature of the crucible is measured by 
means of a thermocouple R located directly beneath 
the crucible, and the output of this thermocouple is 
fed through the potentiometer buckout D to a second 
channel of the recording potentiometer C. Portable 
potentiometer U is used to check the temperature 
of thermocouple R. Thus, a continuous record of 
the erucible temperature, as well as of the weight 
of the sample, is provided on the same sheet. ‘The 
third channel of the recorder is used to monitor the 
temperature of either the lower furnace F or the 
upper furnace E by means of a thermocouple 5, 


The author is indebted to M. Morse of the Electronic Section for his help in 
the design and construction of some of the electronic parts of the apparatus 
‘ Except for some details which are described, furnaces E and F are similar to 
those used in connection with the spring balance in previous work [1]. 
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imbedded in the insulation of E, and connected to 
the recorder through the potentiometer buckout 
D’. Power to the furnaces F and E is supplied from 
the voltage regulators H and H’, respectively, and 
is controlled by the variable autotransformers G 
and G’. The power to the lower furnace is, in 
addition, rectified by the bridge rectifier K. It 
passes through the resistance N, and is controlled 
by the electronic thermostat J to maintain the con- 
stant temperature. The currents fed into F and 
EK are read by means of ammeters M and M’. 

The balance, known as the “Elmic’’ recording 
balance, is a modification of an instrument described 
earlier by Theodore Gast [4]. It makes use of a 
torsion suspension and a feedback circuit with an 
electrodynamometer type of actuation. It has four 
ranges, 10, 5, 2, and 1 mg full seale, the reading 
being indicated on a meter. For continuous re- 
cording, as in the present equipment, a 10-mv 
recorder is connected to the balance. 

The furnace and vacuum system have previously 
been described [1, 2] in detail. The electronic ther- 
mostat [5] consists basically of a resistance-bridge 
sensing unit, the output from which is amplified 
and used with a time-proportional control circuit 
to control the power supplied to the furnace. 

The recorder is a standard 3-channel, 10-mvy, 
commercial type. A photograph of a part of the 
recorder chart is shown in figure 3. The chart gives 
a record of the first 7 hr of operation of an experi- 
ment, which is described later, on the rate of deg- 
radation of polymethylene at 362.4° C. The hori- 
zontal seale in the chart can be used to calculate the 
loss in weight and the temperature of the thermo- 
couples S and R (see fig. 1). The vertical scale shows 
the duration of the experiment. Lines I and I] 
represent records of temperatures indicated by 
thermocouples S and R, respectively. Line III 
represents a record of the loss of weight of the sample 
being heated. In this as well as in all the other rate 
experiments described in this paper, the 5-mg full- 
scale range of the electronic balance was used. 

The procedure used in the present work is in 
general the same as that used previously with the 
spring balance, and has been described in detail 
(1, 2}. In the present method it also required about 
15 min from the time the furnace was raised in 
position, for the thermocouple R to reach the required 
temperature (fig. 3). 


3. Temperature Measurement 


In the previous work on rates of thermal degrada- 
tion of polymers the cooling-end effect of furnace 
F, figure 4, was overcome by adding a smaller 
furnace EK on top of furnace F, and keeping the 
temperature at the interface L, constant by means 
of thermocouples W and W’ in conjunction with 
an electronic thermostat [5] not shown in the figure. 
It was assumed that by this means it would be 
possible to maintain about the same temperature in 
crucible () as that indicated by thermocouple R. 


To check this, a series of blank experiments was 
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carried out in which a chromel-constantan ther- 
mocouple V (shown in fig. 4), was placed inside the 
crucible Q and pressed against the bottom of the 
crucible by means of a platinum ring. The crucible 
was suspended in the glass apparatus in the same 
position as would be used in an actual experiment, 
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Figure 3. Part of recorder chart, indicating the first 7 hr of 


operation of an experiment on the degradation of polyme thyulene 
at 362.4" C. 


Femperature indicated by thermocouple S (1); temperature indicated by 


thermocouple R (11); loss of weight of sample (111 


When operating under the experimental conditions 
previously used, there was a temperature difference 
between thermocouples R and V of about —14° © 
However, it was found that the thermocouples \ 
and R could be made to record the same temperature 
for a distance of 2 to 3 mm between the crucible 
and thermocouple R by adjusting the current input 
in the furnaces E and FF. The system was then 
calibrated by measuring the required currents for 
a series of temperatures and constructing 
responding current-temperature curves in the range 
200° to 400° ©. These curves were followed in the 
experimental work on rates described. Temperature 
control at the interface L (fig. 4) was not used. 


4. Materials Used 


The following materials were used in this investiga- 
tion: 

1. Polwstyrene. Prepared thermally by the Dow 
Chemical Co. Average molecular weight of 230,000, 
as determined by the osmotic pressure method, 
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PiGuri / Details of furnace and control of in 


lemperalure 


Prepared thermally by 


2. Polv(a-methyvistvrene 
low-temperature poly- 


the Dow ( hemical ('o. by 


merization \lolecular welglit by liehit scattering, 
350.000 

Polymethvlene Prepared aecording to the 
method described by Buckley and Ray [6] by 


decomposing diazomethane, using trimethyl borate 


as catalyst Its molecules are unbranched without 
methyl groups, and of high molecular weight [7] 

1. Polvethyvlene 1. A pure grade polymer, highh 
brane hed, and of high molecular welglit Some of 
Its properties are FIVE In table | 
» Polvethvlene 1] A polymer ol high purity, 
having il weight-average molecular welglht of about 
20,000. Its branching characteristics are not known. 

6. Poly chlorotrifluoroetlivlene Prepared with- 
out catalyst and having a weight-average molecular 
weight of 100,000 

All these polymers, except number t, are the 
sume materials as those used previously in work on 


TABLI | Properties of branched polyeth jlene [* 


rates of thermal degradation at higher temperatures 
in the tungsten spring balance apparatus [1,2,7,8,9]. 
As in previous work, the samples were heated in a 
vacuum at 100° C to a constant weight before they 
were used in the rate experiments. 


5. Experimental 


5.1. Rate Measurements 


All the experiments on thermal degradation were 
carried out in a vacuum of about 107* mm of mer- 
cury, obtained by means of an oil pump and mercury 
diffusion pump. Weights of the samples varied from 
ft to 5 mg. The furnace was first heated to the 
operating temperature in preparation for an experi- 
ment, and then raised into position by means of an 
elevator for pyrolysis of the polymer sample. It 
ordinarily took about 15 min from the time the 
furnace was raised into position to the time when 
thermocouple R (figs. 1 or 4) indicated the operating 
temperature. The temperature remained constant 
throughout the experiment to) within 0.1° C 
Zero time for a given experiment was taken as the 
time when this temperature was reached (fig. 3). 
At the end of an experiment the residue in the 
crucible was weighed in order to check on the weight 
loss indicated on the recorder chart. The two values 
differed by about 3 percent. This difference is 
within the combined experimental errors involved in 
weighing small samples, particularly the residues, on 
a semunicrobalance used in this work and in reading 
the recorder chart. The balance weight was used in 
calculating the results and the error was distributed 
over the horizontal scale of the recorder chart 

All these polymers, except number 4, are the same 
materials as those used previously in work on rates 
of thermal degradation at higher temperatures in 
the tungsten spring balance apparatus |1, 2, 7, 8, 9). 
As in previous work, the samples were heated in a 
vacuum at 100° C to a constant weight before they 
were used in the rate experiments 

In figures 5, 6, 7, and 8, percentage volatilization 
is plotted as a function of time for polystyrene, polv- 
(a-methylstvrene), polymethylene, and polvethvlene 
I, respectively. Generally the curves resemble those 
obtained for the same materials in the spring balance 
apparatus, except for the initial parts. Whereas in 
the previous work, which was carried out at much 
higher temperatures, the losses during the heating-up 
period were In some cases as high as 10 percent, the 
curves in figures 5, 6, 7, and 8, show practically no 
such initial losses. 

Curves showing rates of degradation as a function 
of percentage of volatilization for the four polymers 
are shown in figures 9, 10, 11, and 12. The rates 
were calculated from the volatilization-time curves 
and are expressed either as percentege of sample 
(KX;-curves) or as percentage of residue (K»-curves) 
In the case of poly stvrene the K, curve ts given only 
for 299.3° (dash-and-dot line of fig. 9 This curve 
has & maximum at about 30 percent volatilization. 
In the previous work at higher temperatures [1] the 
maxima were between 35- and H-percent volatiliza- 
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TABLE 2.—Rates of thermal degradation of polymers 


Tem- Duration! Total Initial 

Polwmet pera- of ex- volati- rate 

ture periment lized 

( hr Z 7 /min 

200). 7 31.0 23. 3 0. 0090 

> 209. 3 3. 4 77.6 O190 
Polystyrene | 302. 6 27.8 16.4 _ 0230 
SIL.S 29.8 92.1 . 0484 

| 228.8 51.0 27.5 oO110 

Poly (a-methylstyrene 239.4 27.8 37.8 . 0820 
| 247.1 26.0 60. 4 0720 

| 345. 0 55.3 91.5 0080 

' : ; 354.4 ‘ALO 17.6 022 
Polymethyle ) 362.4 37 60.6 0424 
| 367.0 23.8 63.1 . 0094 
350.9 4.5 39.6 O16 
> $55. 4 31.0 ‘SoH . 0260 

vet em - 
Polyethylene I | 303, 2 31.0 4.7 . 0430 
$72. 6 25.0 75.3 0O50 

Polyethylene II 360, 2 31.0 1.8 O81 
Poly (chlorotrifluorethyvlenc $31.8 460 39,2 043 


tion. The K,. curves do not differ much in shape 
from those obtained previously in the spring balance. 
All the pertinent experimental data on rates of 
thermal degradation obtained in the present in- 
vestigation are given in table 2. 

The K, curves for poly(a-methylstyrene) in figure 
10 differ up to about 20- to 25-percent volatilization, 
from the corresponding curves obtained at higher 
temperatures. Here the rates are initially high, then 
drop gradually in convex lines to about 8- to 20- 
percent volatilization, and finally in straight lines in 
the manner of a first-order reaction. In the previous 
work at higher temperatures the rates were initially 
low, then rose gradually, reaching maxima between 
20- and 25-percent volatilization, and finally dropped 
in the same manner as shown in figure 10. It is 
likely that this difference is due not to the difference 
in temperature but to the fact that in the present 
work the rates are much slower and are recorded 
automatically so that the details are more clearly 
delineated and the results are more dependable. 

The K, curves for polymethylene (fig. 11) are 
similar to those obtained previously at higher tem- 
peratures. Polvethvlene I (fig. 12) had not been 
studied previously. It would not be correct to 
compare the present results with those obtained 
previously for polvethvlene II, because the two 
polymers are different in structure. On the other 
hand, a comparison between the K, curves for polv- 
ethvlene IT and for polymethylene shows a striking 
difference between the thermal behavior of these 
two polymers. The polymethylene rate curves 
start at the origin, rise to maxima at about 3- to 
20-percent loss, then fall in straight lines approach- 
ing zero at about 100 percent loss; while the polv- 
ethylene curves start at a fast rate, curve gradually 
in convex form to about 16- to 30-percent loss, and 
then fall in straight lines like the polymethyvlene 
curves. The explanation undoubtedly lies in’ the 
fact that volatilization in both polymers is due to 
fragments breaking off from the ends of the polymer 
chains; and, since the branched polvethylene has 


more ends, it volatilizes faster initially until most of 
the branches disappear. 

The intercepts of the extrapolated straight parts 
of the rate curves K, or K, in figures 9, 10, 11, and 12 
are the apparent initial rates and indicate what the 
actual initial rates would have been were it not for 
rate-disturbing factors, such as rapid initial drop in 
the average molecular weight, presence of foreign 
groups at the ends of chains, or weak points [10] in 
the chains. However, the effects of these disturbing 
factors disappear after degradation proceeds to some 
extent. 

The shape of the polystyrene rate curves (fig. 9) 
indicates that the reaction approaches a zero order. 
The pending benzene rings seem to have a stabilizing 
effect on the chains, so that the number of chain- 
ends and, therefore, the rate of formation of monomer 
dimer, and trimer at these ends by an unzipping 
process stays constant. In the case of poly (a- 
methlystvrene), the presence of quaternary carbons 
in the chain imparts instability to the chain, so that 
formation of new chain ends takes place at random 
and the reaction is of the nature of a first order one. 
In polymethylene and polyethylene the chains also 
break at random and the reaction is first order, 


5.2. Activation Energies 


The activation energies were obtained by plotting 
the logarithm of the apparent initial rates at various 
temperatures for any given polymer as a function of 
the inverse of the absolute temperature. The results 
are shown in figures 13, 14, 15, and 16 for polystyrene, 
poly(a-methylstyrene), polymethylene, and poly- 
ethylene I, respectively. In figures 13, 14, and 15 
are also plotted the results of the previous work 
11, 2, 8} carried out on the same materials at higher 
temperatures. Judging ‘from the position of the 
dashed line (for the previous work) with respect to 
the solid line (for the present work), there seems to 
be a discrepancy in the two sets of temperatures to 
the extent of 17° for polystyrene, 12° for polv(a- 
methylstvrene), and 14° for polymethylene. 

In the present work the temperatures were checked 
by actual measurements, as described earlier and 
illustrated in figure 4. In the previous work it was 
erroneously assumed that there would be a region 
of uniform temperature within 2 to 38 mm distance 
above the thermocouple R (fig. 4). Actually the 
temperature of the crucible above R must have been 
lower than that indicated by the thermocouple by 
the amounts shown above. <A lowering of the tem- 
perature readings in the earlier work will not change 
inany way the shape of the volatilization-time curves, 
rate curves, or the initial rates, but will lower slightly 
the values of the activation energies involved. Thus 
the activation energies of polystvrene, polyv(a- 
methylstvrene), and polymethylene, after correcting 
the temperatures, are 55, 55, and 73 keal/mole, re- 
spectively, as compared with the values of 58, 58, and 
76 keal/mole reported in the previous work for the 
same polymers. The results for branched poly- 
ethylene I, shown in figure 16, indicate an activation 
energy of 64 keal. 
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As a further check on the temperature readings In 
the tungsten spring balance, two polymers, polv- 
chlorotrifluoroethylene) of average molecular weight 
weight 100,000, and polvethyvlene IL of molecular 
weight 20,000, were pyvrolyzed in the electronic bal- 
ance at relatively low temperatures. The results 
are shown in figures 17 and 18. Pertinent data for 
these polymers are given in table 2. The pyrolysis 
curves and the rate curves resemble closely those 
obtained previously |1, 9| for these two polymers at 
higher temperatures. However, when the results 
obtained in figures 17 and 18 and those obtained in 
the previous work are plotted in terms of logarithm 
of the initial rate as a function of inverse of absolute 
temperature, we find that the temperature readings 
in the previous work are too high by 14° and 13° C 
for poly (chlorotrifluoroethy lene and poly ethy lene | |, 
respectively, and the activation energies too high by 
3 keal/ mole for both polymers. 


5.4. Correction of Previous Work 


500 385 320 TEMPERATURE On the basis of data in table 2 and the time and 
= v rate curves shown, it is suggested that corrections 

5 © ¢ be made in the temperature data and activation 
energies obtained by the present author and his CO- 
workers in their studies of rates of thermal degrada- 
tion of a number of poly mers determined prey ioushy 
‘ in the tungsten spring balance. The suggested cor- 
rections for temperatures and the corrected activa- 

tion energies are shown in table 3. The corrections - 

given in this table for polystyrene of molecular weight 

230,000, poly (a-methylstvrene), polvmethylene, poly- 

. (chlorotrifluoroethvlene . and polvethvlene of mo- 
lecular weight 20,000, are those found experimentally 
in the present Investigation The corrections for 
the other polymers shown in the table are those 
58 based on the average decrease of temperature and 

; activation energies found experimentaliy in this work. 
These averages are 14° C and —3 keal/mole. No 
g corrections are required for Teflon |10| since the 
polymethylene. temperatures for this polymer were calibrated by 

Dashed line indicates previous work comparing the results obtained in the spring balance 
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1 percent of sample per minute 


with those obtained by a pressure method. In the 
more recent work in the spring balance {11, 12] the 
furnace was calibrated in the same manner as de- 
seribed in the present work. 


6. Comparison of Rates Obtained by 
Several Investigators 


In view of the difficulties involved in measuring 
temperatures in a study of thermal degradation of 
polymers, there is a likelihood that there should be 
disagreements in the results obtained by various in- 
vestigators. In figure 19 the results of rate studies 
obtained by Jellinek [13,14], Grassie and Kerr [15], 
and the present author are compared. The results 
of all these investigators were obtained in vacuum, 
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PABLE 3 Corrections to be applied to experimental te m perature 


and to activation energies in previous work 


Experi- 
remper mental and 
Polymer Reference iture corrected 


correction wtivation 


enereles 


( 1 mole 

Polystyrene (230,000 [ 17 
Polystyrene (106,000 | 14 
Polystyrene (584,000 {1 14 
Polyethylene (20,000 | 13 63 
Poly (m-methylstyrene {2 14 oh 
Poly (a-deuterost yrene [2 14 
Poly (a-methylstyrene [2 12 
Hivdrogenated polystyrene {2 14 1 
Poly (methyl methacrylate) A {2 14 i) 
Poly(methyl methacrylate) B [2 14 2 
Poly (methyl acrylate {2 14 $4 
Polymethylen [s 14 2 
Polypropylene |S 14 s 
Polyisobut viene [s 14 10 
Polybenzy| {sS 4 0) 
Poly (8-deuterostvrene high mole 

cular weight [s 14 Mw 
Poly (chlorotrifluoroet hy lene ’ 14 7® 
Poly (a@,8,8-trifluorost yrene 1 4 4 
Poly (p-xyly lene {y 14 73 
* In reference [9] the tivation energy for poly(chlorotrifluorcethylene) is 


given as 66 keal/mole his is a misprint and should have read 60 Kea When 


educed by 3 units the activation energy becomes 


using polymers prepared without catalysts and hav- 
ing high molecular weights. The solid lines in figure 
19 represent polystyrene, and the dash-and-dot lines 
represent poly(a-methylstyrene); also, the letters J, 
G, and M, in front of the temperatures indicate the 
work of Jellinek, Grassie, and Madorsky, respec- 
tively. Temperatures in parentheses represent the 
uncorrected temperatures obtained by the present 
author in the tungsten spring balance, 

The polystyrene curves M-290.7° and M-311.8 
are the same as those given in figure 5. Curve 
G-290° was obtained by Grassie and Kerr [15] on 
a sample of polystyrene (S—3), prepared without a 
catalyst and having a molecular weight of over 
300,000. Their rate apparatus consisted of a dy- 
namic molecular still heated electrically [16]. The 
G-290° curve is very close to the M-—311.8° curve, 
so that the temperature given by Grassie and Kerr 
is at least 20° lower than the one obtained in the 
electronic balance in the present work. The \1-290.7 
curve lies far below the G-—290° curve. At higher 
temperatures the J-364° curve was obtained by 
Jellinek for a polystyrene sample (Fir), prepared 
without catalysts and having a molecular weight of 
200,000 {13, 14]. He used a quartz spring balance. 
Curves M-328°(345°) and M-848°(365°) were ob- 
tained by the present author on a thermally prepared 
polystvrene of molecular weight 230,000, using a 
tungsten spring balance [1]. Position of J-364° is 
close to that of M-—328°(345°) and is far below that of 
\M-348°(365°) curve. Thus Jellinek’s temperature 
reading is about 30° to 32° higher than that of the 
present author and about 50° higher than that of 
Grassie and Kerr. 

In the case of polv(a-methylstvrene ) Jellinek {1 1 | 
pyrolyzed a 70,000 molecular weight sample in the 
quartz spring apparatus and obtained a volatiliza- 
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Figure 10 Comparison of results obtained by various inves- 
tigators on thermal degradation of polymers at various ltem- 
peratures 
Polystyrene poly (a-methylist vrene Che letters 
in front of each temperature refer to the investigators: G, Grassie: J, Jellinek 
M, Muadorsky remperatures in parentheses are uncorrected temperatures 
obtained the spring-balance pyrolysis apparatus 
tion-time curve shown in figure 1S as J-28] The 
curve M-247.1° is the same as the 247.1° curve given 
above in figure 6. Jellinek’s temperature reading 


appears to be about 34 higher than the one obtained 
by the present author using the electronic balance. 

It is believed that the principal reason for the 
differences between the various investigators lies in 
measurement of the temperature of the sample. As 
shown earlier in this paper, this measurement is very 
sensitive to the maintenance of a uniform te mpera- 
ture in the zone where the sample and the thermo- 
couple are located and the relative position of 
these two in this zone. 


to 


7. Discussion and Conclusions 


A quantitative determination of rates, and conse- 
quently of activation In the study of 
thermal degradation of polymers by the pressure or 
loss-of-weight method i by many difficulties 
The following are some of the more serious difficulties. 

1. Spattering caused by the bursting of bubbles of 
gaseous products of degradation which accumulate in 
the polymer sample 

2. Delay in the vaporization of 
due to slow diffusion through the 
polymer. 


energies, 


s beset 


volatile products 
solid or semisolid 
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The difficulty of maintaining a uniform temper- 
ature through the polymer sample due to its poor 
heat conductivity. This is further aggravated by the 
fact that the degradation reactions involved are 
endothermic in most 

The difficulties associated with measurement of 
temperature. 

5. Uncertainties in temperature and rate observa- 
tion during the initial heating-up period. 

With regard to points 1, 2, and 3, the difficulties 
were reduced to a minimum in the present work as 
well as in the earlier work by limiting the weight of 
the samples to about 5 mg or As to point 4, 
precautions were taken in the present work to obtain 
a uniform temperature in the region where the 
thermocouple and the sample are located in close 
proximity to each other. The difficulties encountered 
during the initial stage of degradation, referred to in 
point 5, were reduced in this work to a minimum by 
operating at low temperatures, so that very little loss 
bv volatilization took place during the 15-min heating- 
up time 

A study of rates of degradation at low tempera- 
tures, in conjunction with automatic recording of 
loss of weight and of temperature, not only renders 
the. observations more accurate, but also reveals 
details that are otherwise masked by a rapid rate of 
volatilization at high temperatures. This is partice- 
ularly true for polyethylene and polv(a-methyl- 
stvrene) where the actual initial rates of degrada- 
tion are very high. 


Causes. 


less, 
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P s in Metals Caused by Collision With Liquid Drops 
and Soft Metal Spheres 
Olive G. Engel 


An equation is developed to give pit depth as a function of collision velocity for pits 
formed in soft to medium-hard metal plates as a result of collision with liquid drops. The 
rear face of the target plate must be a free surface. The plate thickness must not be less 
than 1.5 to 2.0 times the drop diameter nor greater than 4 to 5 times the drop diameter. _ It 
is shown that, under the same conditions on the target plate, the equation is also applicable 
to pits formed in soft to medium hard metal plates as a result of collision with spheres of 


the same metal that flow radially as a result of the collision. 

data obtained in other laboratories were used to test the equation. 
were copper, 1100—0 aluminum, 2024-0 aluminum, lead, steel, soft iron, and zine. 
was used as the drop liquid against copper, aluminum, lead, and steel. 


Pit-depth-versus-velocity 

Metals used as targets 
Mereury 
Water was used 


as the drop liquid against copper, aluminum, and lead. Spheres of copper, aluminum, 
lead, soft iron, and zine were used against targets of the same materials, respectively. The 
equation can be used to calculate the dynamic compressive vield strength of soft to medium- 


hard metals. 


1. Introduction 


Current research on pits produced in high-speed 
collisions includes: (a) the investigation of solid 
projectiles impinging against solid targets (artillery 
experiments); (b) the investigation of solid projectiles 
impinging against liquid surfaces (water entry 
problems); (c) the investigation of liquid drops 
impinging against solid surfaces (high-speed rain- 
erosion research); and (d) the investigation of liquid 
drops colliding with liquids. The work that is 
described in this paper was initiated as part of a 
high-speed rain-erosion research program and_ the 
entire program was conducted under the sponsorship 
of the Materials Laboratory, Directorate of Labora- 
tories, Wright Air Development Center. 

It is difficult to test for the rain-erosion resistance 
of structural materials at verv high collision velocities 
because of the problems involved in accelerating 
either waterdrops or test specimens to the velocities 
in question. It has been suggested that it may be 
possible to bypass these problems by using drops of 
high-density liquids instead of drops of water for 
the erosion tests.' To develop this idea into a 
reliable test procedure, it is necessary to know the 
corresponding velocities-for-equal-damage when the 
test specimen collides with a drop of a high-density 
liquid and when it collides with a drop of water. 
The pit-depth-versus-velocity equation presented 
in this paper was developed to provide this 
information. 

It has also been found that this pit-depth-versus- 
velocity equation is applicable without change to 
collisions of spheres of the soft, ductile metals with 
targets of the same metal in those cases in which 
the sphere flows as a result of collision with the metal 
plate when the collision velocity is as high as 5,000 
Pellets of the soft, ductile metals appear, 


This idea was suggested by the rain-erosion research group working at Convair 
Division of General Dynamics Corp. in San Diego, Calif 


ft see. 
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therefore, to behave as though they were liquid drops 
when they impinge against a solid surface at this 
velocity. 

The experimental work reported in this paper was 
done in other laboratories.” 


2. Collisions With Liquid Drops 


The damage done to solid materials that collide 
at high speed with liquid drops depends both on 
the properties of a drop of liquid in collision and on 
the characteristic properties of the solid. 


2.1. Damaging Properties of a Liquid Drop in High- 
Speed Collisions 


As the relative collision velocity changes, the 
effects of collisions with a liquid drop also change. 
In high-speed collisions with the planar surface of a 
solid the liquid drop acts as though it were a hard 
solid sphere, but, unlike a sphere of hard solid 
material, it undergoes an ultrarapid radial flow 
outward about the point of impingement [1,2]. 
These damaging properties of an impinging liquid 
drop vary in intensity depending on the density of 
the liquid of the drop, on the relative collision 
velocity, and on the extent to which the solid surface 
vields under the blow. At any given impingement 
velocity, collision of a ductile metal plate with a 
liquid drop does not produce as deep a pit as collision 
of the same metal plate with a hard solid sphere 
would produce. This is because when a planar 
solid collides with a liquid drop part of the collision 
energy is transformed into the radial flow of the 
liquid. 


2 The pit-dlepth-versus-velocity data reported for collisions of solid targets 


with liquid drops were obtained by the rain-erosion research group working at 
Convair Division of General Dynamics Corp 
The pit-depth-versus-velocity data for collisions of spheres of the soft metals 
that flow as a result of the collision were obtained by Partridge, VanFleet, and 
Whited at the University of Utah (see reference [11 
‘ Figures in brackets indicate the literature references at the end of this paper. 








2.2. Response of Solid Target Materials 


The extent and type of damage produced on solids 
as a result of collision with liquid drops depends 
strongly on the characteristic properties of the 
solids [1]. If the damage is to be described mathe- 
matically, it will be necessary to develop separate 
equations for the damage produced on solid mate- 
rials that have widely different properties. 


2.3. Dimensional Analysis of the Damage 


To develop an equation that will give the depth 
of pits produced in high-speed collision between the 
planar surface of a solid and a liquid drop is a 
difficult) problem from a theoretical standpoint. 
Dimensional analysis has proved to be a_ useful 
tool in the solution of difficult problems in the past. 
It is used in this paper to develop il pit-depth- 
versus-velocity equation. 

The method of dimensional analysis has been 
discussed by Buckingham [3, 4], Bridgman [5], 
Birkhoff iG], and others. If a phir sical process Can 
be deseribed by phvsical quantities of n different 
kinds, and if none of the quantities involved in the 
has been then the process 1s 
deseribed by the equation: 


process overlooked. 


Pirin Wiles cs x Wes Pty Vey 8 es (), 
where the (’s are the physical quantities involved 
in the process, and the r’s are ratios. The ()’s are 


quantities of different) kinds. If several quantities 
of anv one kind are involved in the process, they are 
specified by the value of any them, and by 
ratios of the others to this one. These ratios are 


one of 


the 7's of eq (1 
A certain number & of the @’s are 
fundamental. The & selected @’s comprise a possible 
fundamental dimensions and the remaining 
(/-quantities can be expressed in terms of them 
After the /& fundamental @’s, now designated as 
S’s, have selected, the remaining V's are 
designated as P-quantities. In a mechanical system 
the number of the S-quantities is three because the 
total number of dimensions required to express any 
of the Q-quantities involved in a mechanical process 
is three. These dimensions are usually mass, length, 
and time 
Kquation (1 


selected as 


set of 


been 


may be put in the form 


F(a, te, o~«e ® aS Oe Re (), 2 


where the w’s are independent produets of the argu- 
ments (Y, and are dimensionless in the fundamental 
units, or in the form 


in which aw, is any of the z’s. 
A x-product is formed by multiplying each /’- 
quantity by the S-quantities, which are 


raised to 


powers a, 8, x are required so that 
cancel those of the P-quantity 
That ts, 


whatever 
their units will just 
to make the z-product dimensionless. 


a 5 « > 
rT (Qt-Qe-... ODP. (4) 


From a consideration of the equation for the 
pressure that results when a solid surface runs into 
a liquid drop at high speed [2], of the expression for 
the radius of flow of the liquid about the central 
point of the collision [2], and of the well known 
equation’: for the shear between lavers of 
liquid in laminar flow, the quantities (see table 1 
that should prove to be important in determining 
the damage that results from the collision of a solid 
surface with a liquid drop are: ¢, ¢’, p, p’, Vi wu, ¥, Y'; 
d, and 6’ where ¢ is the speed of sound, p is the dens- 
itv, Vis the relative impingement velocity, u is the 

viscosity, yis the surface tension, d is the diameter of 
the liquid ‘drop, and 6’ is the damage parameter. 
Primed qui intities refer to the solid material; 
unprimed quantities refer to the liquid of the drop. 
In addition to these quantities there is the energy 


st ress 


per unit volume put in, /’, and the energy per unit 
volume returned, e’, during the collision. £7 is the 


energy per unit volume that the solid material can 
absorb without nonrecoverable deformation or frac- 
ture; ¢’ is the energy per unit volume that the solid 
material can return. The ratio of these quantities, 
e’/I-’, is a measure of the resilience of the solid 
material. 


TABLE l. The three fundamental S-quantities, the P-quantilties, 
and the r-ralios involved mn ordain erosion damaude 
) I) ! S ] Dimer 
| et | L/7 = 
Drop d l r / S 
1) ty iquid Vi S 
\ t | VL / i\ 
Re l 
Numt 
~ t V7 / pid \ 
juid Wet 
Nu t 
\ i Lo7 / \ NI 
qu N } 
liq 
Da I ! l / 
Yield or ru i VL i Tr I \ 
ergy det 
lid. F 
K 1 
I 
R: wrfev 
| I 


Of these quantities V, d, and p are arbitrarily 


selected to be the S-quantities, S,, Sj, and Ss, 
respectively. Considering both these S-quantities 
and the remaining quantities, there are four pairs 


of quantities of the same kind: ¢, ¢’; p, as 4% 
and e’, &’. These are expressed as four 
ri, Yo, 73, and ry, respectively. One member of each 


of these pairs is a ?-quantity if it has not already 


ratios, 
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been chosen be an S-quantity. Also any other 
of the tabulated quantities that has neither been | 
expressed } in a ratio nor chosen to be an S- quantity 
isa P?-quantity. The ?-quantities are then: 
u, y, c, 6’, and E” (see table 

The z-products are formed by multiplying the 
P-quantities by the product of powers of the S- 
quantities. Hach S-quantity is raised to whatever 
power is required to make the z-product dimension- 
For example, the viscosity uw has the dimen- 
M/(LT) where AM is mass, L is length, and 7 


is time, and, therefore: 


less. 


SLOnS 


(M/Z) (L) (L/T) (LT/M). (5) 


rw, pd Vip 


Surface tension y has the dimensions 7/7? and, 


therefore, 

wo '=pd V?/y=(M/L*) (L) (£7/T*) (T?/M). (6 
m;, 7, and ws, are formed by a similar procedure. 
mr, ‘and wr‘ are the dimensionless Reynolds Number 
and the square of the Weber Number, respectively. 
mr; 'isakind of Mach Number that gives the ratio 
of the relative impact velocity between the — | 
drop and the solid surface to the speed of sound i 
the liquid of which the drop is composed. 

The damage parameter 6’ may be chosen to be a | 
length, a volume, or an area. A length is the easiest 
quantity to measure experimentally. Taking 6’ to 
be a length, for the damage caused by high- 
speed collision of a solid surface with a liquid drop is 





eq (5 


5’ /d—F'[pdV'/u, pdV2/y, Vie, E’/(pV2), e’/e, 


Kk’). (7) 


p/p, v'/y, ¢ 


It was pointed out in section 2.2 that it will be 
necessary to develop separate equations to describe 
the damage produced in liquid-drop-solid-surface 
collisions for materials that have widely different 
properties. The damage produced on nonrubbery 
materials that fail by plastic vield but without frae- 
ture the to consider. It ineludes 
materials such as the soft and medium hard metals, 
If the collision velocity is not too high, the damage 
marks produced on these materials by high-speed 
collision with liquid drops are empty spheroidal seg- 
ments. Surface and cross-sectional views of such 
pits in collision experiments performed elsewhere ? 
are shown in figures 1, 2, and 3. 

The damage pits produced in the planar surface of 
soft and medium hard metals when the leading sur- 
face of the solid runs into a liquid drop at high speed 
are almost exclusively due to the impact pressure 
produced. Although the radial flow of the liquid 
does make a contribution to the damage in the case 
of the soft metals (it drags metal up the walls of the 
damage pit and piles it up at the mouth of the crater), 
this contribution may be small enough in comparison 
with the damage that is caused by the pressure so 


Is easiest case 


that it can be neglected in first approximation. 
Neglecting all the dimensionless products that result 
from consideration of the radial flow of the liquid 
drop and from resilience of the solid target material, 
the approximate dimensionless damage equation for 
this class of materials is 

5’ /d kh’ 


F'[p’/p, e’/e, Vie, (p\"?)} 


where 6’ is the depth of the damage pits. 
2.4. Experimental Observations 


values of pit depth plotted 
against relative collision velocity indicate that in 
the velocity range that has been investigated the 
pit-depth-versus-velocity curve is very close to a 
straight line regardless of whether the pits were 
caused by firing steel spheres against a stationary 
specimen plate or by firing a specimen plate at high 
velocity against re latively stationary liquid drops. 
It is to be expected that at a fixed relative collision 
veloci ‘itv the de pth of pit caused by a a projec ‘tile that 
does not flow during the collision (hardened steel 
sphere) will be deeper than that caused by a projec- 
tile that does flow during the collision (liquid drop). 

In the first case all of the kinetic energy (neglecting 
that which is converted into heat) is delivered to the 
solid target; in the second case part of the kinetic 
energy is used to produce the flow of the projectile 
and only the remainder is delivered to the solid 
target. 

Experiments in which 7/32-in. steel spheres were 
fired against 1/8-in.-thick type-1100 aluminum 
plates have been carried out by Mr. Herschel L. 
Smith at the Naval Research Laboratory, Washing- 
ton, D.C. Some of the data obtained in these ex- 
periments are presented in figure 4. One important 
fact is apparent from these data, namely, for the 
sphere diameter and for the thickness of specimen 
plate used, if the specimen plate has only peripheral 
support the pit-depth-versus-velocity curve has 
an intercept on the velocity axis, but if the specimen 
plate is rigidly backed by a heavy steel supporting 
plate the pit- -depth-versus-velocity curve goes directly 
to the origin. 

Experiments in which target plates of different 
metals were fired into drops of mereury and drops 
of water at high speed have been carried out in an- 


Graphs of measured 


other laboratorv.2. Some of the data obtained are 
plotted figures 5 and 6. It can be seen that in 
every case the pit-depth-versus-velocity curve. is 


essentially a straight line with an intercept on the 
velocity axis. From the bulge on the trailing face 
of the specimen in picture 6 of figure 2 it is apparent 
that the metal plates are supported during the 
collision with the liquid drops in such a way that 
the trailing face of the target plate is a free surface. 
An important fact apparent from the data plotted 
in figure 6 is that the slope of the straight line is 

function of the drop size. 

On the basis of the evidence presented, if the speci- 
men has peripheral support only, the pit-depth- 
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SURFACE VIEW 4 CROSS—SECTIONAL VIEW 
IMPINGEMENT VELOCITY 686 ft/sec 





2 SURFACE VIEW 5 CROSS—SECTIONAL VIEW 
IMPINGEMENT VELOCITY 1250 ft/sec 





4% SURFACE VIEW r CROSS—SECTIONAL VIEW 
IMPINGEMENT VELOCITY 1645 ft/sec 


kicurReE 1 Vi lrop damage pits in lead, 
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SURFACE VIEW 4 CROSS— SECTIONAL VIEW 
IMPINGEMENT VELOCITY 686 ft/sec 





2 SURFACE VIEW 5 CROSS— SECTIONAL VIEW 
IMPINGEMENT VELOCITY 1285 ft/sec 








3 SURFACE VIEW 6 CROSS—SECTIONAL VIEW 
IMPINGEMENT VELOCITY 2320 ft/sec 


FIGURE 2 Vercury-drop damage pits in 1100 aluminum, 
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SURFACE VIEW 4 CROSS—SECTIONAL VIEW 
IMPINGEMENT VELOCITY 695 ft/sec 


2 SURFACE VIEW CROSS —SECTIONAL VIEW 
IMPINGEMENT VELOCITY seo ft/sec 











3 SURFACE VIEW 6 CROSS—SECTIONAL VIEW 
IMPINGEMENT VELOCITY 2445 ft/sec 


FiGcurRE 3 Vere j-drop damage pils in copp 
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versus-velocity equation for materials of the response 
t\ pe beme considered appears to be: 
i’ =K,V—K,, 9 


where Ay and A, are constants.’ It appears that 
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FIGURE 6. 


0.10-cm mercury drops; A 0.20-cm mercury drops 


this equation applies regardless of whether the pits 
are produced by the impingement of solid projectiles 
against solid targets or by firing solid target plates 
into liquid drops at very high velocities. However, 
because the amount of the kinetic energy delivered 
to the target depends upon whether or not the pro- 
jectile flows during the collision, the constants 
A, and A, will be different for solid-target-versus- 
solid-projectile and for solid-target-versus-liquid- 
projectile collisions. 


2.5. Pit-Depth-Versus-Velocity Equation for Medium- 
Thin Metal Targets 


For targets that are metal plates, that have a 
thickness several times the diameter of the projectile, 
and that are mounted so that the reverse side of the 
plate is a free surface, a pit-depth-versus-velocity 
equation can be developed by considering the move- 
ment of a cylindrical core of material through the 
target plate under the area of contact involved during 
the collision. If the target plate is fired against the 
projectile, this core of target material is slowed down 
with respect to the remainder of the target plate as 
a result of the collision and there is a relative motion 
between this core and the remainder of the target 
plate, which moves forward with respect to it. 


Partridece, VanFleet, and Whited [11] also reported a linear relationship be 
tween pit depth and velocity; the linear relationship was established independ 
ently by the author on the basis of the Convair data 
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Jf the reverse side of the target plate is a free 
surface and the core can move freely with respect 
to the remainder of the plate, the kinetic energy 
transformed as a result of the collision is largely 
converted into the work of moving the core against 
the bonding that holds it to the remainder of the 
target plate lor conditions there is) an 
intercept on the velocity axis. In terms of this model, 
it would appear that the intercept velocity is related 
to the shear vield strength of the material of the 
target plate. In the range of velocities that are just 
in excess of the threshold velocity required to produce 
any permanent damage at all (intercept velocity), it 
appears that the core may only be displaced with 
respect to the remainder of the target plate. For 
this case plastic deformation occurs only at the 
evlindrical boundary of the core. At higher velocities 
plastic deformation may occur within the core itself. 

Consider that a small target plate moving at 
velocity Vin a stationary coordinate system, which 
| target plate, strikes a 

The result of the 


these 


is located outside of the 
stationary liquid drop (see fig. 7). 
collision is that a wave of compression ts initiated 
both in the solid material of the target plate and in 
the liquid of the drop It is helpful to view the 
collision incident as it is seen by an observer located 
at the origin of the coordinate system in figure 7 
This observer sees a target plate, all parts of which 
are moving forward uniformly at velocity Vin the 

y)-direction, approach a spherical drop, all parts 
of which are stationary After the collision 
occurred, the observer zone in which the 
particles of the material of the target plate have 
y)-direction where 


has 


sees a 


taken on a velocity, aV, in the | 
aisa constant; this zone is indicated schematically 
with a dotted line and the letter A in figure 7. 
Similarly, the observer sees a zone in the liquid drop 
in which the particles of the liquid have taken on a 
velocity, 4V. in the (+ y)-direction where 6 ts 
a different constant; this zone is indicated schemati- 
cally with a dotted line and the letter B in figure 7 
Zone A spreads at the characteristic sound velocity 
of the target material, ce’, through the thickness, / , of 
the target plate and zone B spreads at the charae- 
teristic sound velocity, ce. of the licyuid of which the 
diameter, d, of 


drop is composed, through the 
the liquid drop 

To the observer who looks first at a pomt 1 that 
is in the target plate but outside the spreading 
boundary of zone A and who then looks at the mate- 
rial within zone A, it appears that the particles at 
A’ have a velocity Vand that the particles in zone A 
have a velocity (l-a)V. After a time ¢ that is just 
long enough for zone A to complete one trip through 
the thickness of the target plate, the leading surface 
of the target plate at point .1’ in figure 7 has ad- 
vanced a distance Vt, but in the same time f¢ the 
area of the leading surface of the target plate that 
constitutes the upper boundary of zone A in figure 
7, has advanced only by (l-a)Vt. If the average 
negative velocity of the core is taken to be aV)2, 
the displacement that exists between the leading 
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FIGURE 7 Stationary coordinate system for a liquid-drop 


ver sus-solid-specime n collision 


surface of the target plate at point A’ and the upper 
boundary of zone A after time ¢ is aV7/2. 

If the observer watches zone B as it just spreads 
to the opposite side of the liquid drop, he sees that 
this zone reflects from the free air-liquid interface 
with change of sign, that is, he sees a tension wave 
begin to move toward him in the (—y)-direetion 
of figure 7. Within the tension wave the particle 
velocity is 2AV' in the ry -direction. The tension 
wave moves toward the observer at the characteristic 
speed of sound, c, of the liquid of which the drop is 
composed. In the zone of the liquid of the drop 
that has been traversed I the tension wave the 
particles of the liquid are unstressed and are mov Ing 
at velocity 26° in’ the y)-direction. At the 
Instant that the tension wave returns to the collision 
surface, where radial flow of the liquid of the drop 
is occurring, all of the liquid of the drop is unstressed 
and moving at the velocity 26° in the (+ y)- 
direction. 

It can be assumed that return to the collision 
surface of the tensile wave in the liquid provides a 
cutoff for the collision. At this instant all of the 
liquid of the drop is unstressed. The collision sur- 
face will appear to be a free surface to the com- 
pressed particles of the material of the target plate 
in zone A. A wave of relief (tension) will then be 
initiated in the material of the target plate and will 
begin to trail zone A which continues to spread 
through the thickness, d’, of the target plate 
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Let the arbitrary time ¢ referred to previously 
now be the time that is required for zone B to make 
one round trip through the liquid drop. Then 
t—2d/e. The distance to which zone A has spread 
in this time is c’t and this determines the length of 
the compressed zone A which continues to move 
through the target plate. When zone A reaches the 
free trailing surface of the target plate, it also 
reflects as a tension wave. This tension wave moves 
in the (-+-y)-direction of figure 7. In that part of 
zone A that has been traversed by the returning 
tension wave the particles of the material of the 
target plate are under zero stress and have velocity 
(l1—2a)Vo in the (+y)-direction. When the re- 
flected tension wave in the target plate reaches the 
leading surface of the plate it undergoes partial 
reflection back into the material of the target plate 
and partial transmission into the liquid of the drop 
which is running off radially around the central point 
of the collision. 


a. Slope of the Pit-Depth-Versus-Velocity Curve 


The compressional wave spreads slightly as it 
moves through the thickness of the target plate; 
consequently, the core of material through the target 
plate that is slowed down as a result of the collision 
is in reality somewhat conical in shape. For sim- 
plicity, the true situation may be idealized in two 
wavs. First, the may be regarded as a true 
cylinder which is free to move or slide in the diree- 
tion of the collision blow but which is restrained 
laterally. A similar evlinder exists in the liquid of 
which the drop is composed. Secondly, the compres- 
sional waves that move through the eylinder in the | 
target plate and through the evlinder in the drop 
may be regarded as plane waves. 

For the of plane waves, a=2/(2 and 
b { ‘ where Cp and = 2” c’ p’ are the 
acoustic impedances of the material of the drop and 
of the material of the target plate, respectively. 
The relations are found by equating the plane-wave 
stress in the material of the drop, ¢, to the plane-wave 
stress in the material of the target plate, 7, that is, 


core 


Cruse 


o—cpr—c' pt w; (10) 
where ¢, 7’ are the particle velocities in the zones 
traversed by the compressional waves initiated in 
the drop and in the target plate, respective ly, and 
bv imposing the condition that the impacting surfaces 
remain in contact, namely, 


o+e’=—V. (11) 
Krom eq (10) and (11 
py’ cov c’p’) 
-V/( 4 
aV 
and o=2 V/( ) 
hV (12) | 


\ 


This derivation is the same as that for the collision 
of two free rods except that, for the case that the 
rods are cylinders that are free to move in the colli- 
sion direction but that are restrained laterally, the 
sound speed that must be used is the speed of sound 
be 


in an infinite medium. This can shown as 
follows. According to Hooke’s law: 
Ke, =o; —v(20,) =o; —2va,, (13 


where / is Young’s modulus, e; is the longitudinal 


strain, o; is the longitudinal stress, v is Poisson's 
ratio, and is the radial stress. Also, 
Ive, Ce" WG2 “TF Gz) 0, (14) 
where ¢, is the radial strain. From eq (14) 
o,=—vo,/(l—vp). (15 
By substituting the expression for o, given by eq 
(15) into eq (13) it is found that 
(l—v) Ke, 
07; ? 16 
(1—2y)(1+-p) 
and because, for plane waves, c=|o,/pe,|'” 
(l—yp) |!” 17 
(l—2y)(1+y) p 


This is the sound speed for an infinite medium and 
it is the sound speed that must be used both for ¢ 
and tor ¢c’ in the expressions for a and for 6. 

The depth of pit produced at any relative impinge- 
ment velocity is proportional to the negative velocity 
produced in the evlindrical core in the target plate 
as a result of the collision and to the time that this 
negative velocity exists. If it is assumed that the 
elastic wave that is induced in the evlindrical core 
asa result of the collision makes a round trip through 
the core in the time interval 2d/¢, then, at the end 
of this time interval, the cylindrical core is moving 
as a rigid body and the negative velocity that it has 
acquired with respect to the remainder of the target 


plate is given approximately by 22V7/(2+ 2’). During 
this time interval the average negative velocity given 
tothe material of the evlindrical core was 2V/(2+ 2’). 


The slope A, of (9) is then &,(2d/c) [2/(2+2’)} 
where A, is a new constant, 

The experimental data for a number of 
metals of widely different: properties fired 
drops of mercury, 


eq 


target 
against 
and for three target metals fired 
both against drops of mercury and drops of water, 
require that k;—3.6. Beeause this value of the 
constant ky brings the calculated curve into good 
agreement with the experimental data for a number 
of target metals and for two drop liquids, it would 
appear that &,; does not involve properties either of 
6 The author is indebted to Dr, John M,. Frankland of NBS Mechanics Section 


for the proof in elasticity theory that it is the sound speed in infinite medium that 
is required 
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the material of the target plate or of the liquid of the 
drop. Therefore, eq (9 is 


: bul, 18) 


; 


Substitution of - Cp, z cp 
eq (18) in the dimensionless form 


‘and division by d puts 


«/ - 
‘ 


to 


p 
where A,/d is a dimensionless intercept. Four of the 
dimensionless quotients predicted by eq Ss)jare found 
in eq (19 


b. Intercept on the Velocity Axis 


The process previously postulated, namely, that 
the core of target material below the collision area 
is slowed down with respect to the remainder of the 
target plate so that the remainder of the target plate 
moves forward with respect to it, certainly occurs 
at all impingement velocities. However, no 
servable damage pit is produced at all below a 
critical velocity, and, if the liquid of the drop Is 
kept constant, this critical velocity is different for 
each of the metals for which data are available. 
(See fig. 5 The intercept velocity is also different 
on the same metal if the liquid of the drop is changed. 

The fact that no permanent damage pit is formed 
below a given velocity which is characteristic for 
each metal, indicates that for velocities lower than 
this critical velocity the relative motion between 
the core and the remainder of the target plate Is 
completely elastic and that no permanent shear 
deformation occurs in the material around the 
evlindrical wall of the core. The intercept of the 
pit-depth-versus-velocity curve appears to be a 
function of the shear vield strength of the material, 
that is, of the remaining dimensionless quotient 
Ie’ /(pV*) of eq (8) taking /’, the energy per unit 
volume absorbed by the target material without 
fracture or plastic vield, to be the shear vield 
strength. 

It has been pointed out that the experimental 
pit-lepth-versus-velocity curve is) a_ straight-line 
function of the velocity. Therefore, the dimension- 
less quotient /” (pV?) cannot be used tor the inter- 
cept without climinating the factor 1/V?. This ean 
be accomplished by multiplving the dimensionless 
quotient ky’ (p\? by the square of the dimensionless 
quotient lve of eq (S) to obtain fe’ ( pc* If this 
expression is substituted for the dimensionless inter- 
cept A, /d in eq (19) and if the Intercept condition 
that 6°/d¢d=0 is imposed, the expression that is found 
for the intercept velocity is not able to account for 
the observed ratio of experimental intercepts. [Lt was 
found that the observed intercept velocities can be 
accounted for if the dimensionless quotient /” (pc? 
is multiplied by the dimensionless quotients 
' ‘, and by a dimensionless numerical 


; ] 
(pp a , Cr 


constant, k., having the value 136.8. Substitution 
of the resulting dimensionless quotient, 136.8 
| age (p' 29" Acc’ 2’), for the dimensionless intercept, 
K,/d, in eq (19) and use of the acoustic impedance 
for ep produces the dimensionless pit-depth equa- 
tion 
1I36.8E" 2 
ee : 20) 


‘i pe 


p ‘“p “CC 


If the intercept condition that 6° d=—0 is imposed, 
the intercept velocity, V.. is found to be civen by 


, 19E’(z+2’) 
| . yar 1/9 - (21) 


(pce < ” 


It is noteworthy that one would expect \’, to be 
given by E’ (2+-2’ ‘multiplied by some constant, 
because the plane wave stress is given by 22’V 
(2+-2’). It was found, however, that this expected 
expression would not give the observed intercept 
velocities for both mercury drops and waterdrops if 
the same numerical constant was used for both drop 
liquids. Equation (21) appears at present to be the 
best expression for the intercept velocitv. When 
more data are obtained and the problem is studied 
further, it may be found necessary to change it. 

The numerical constants in eqs (20) and (21) have 
been chosen to vive best fit to the pit-depth-versus 
velocity data. The values used for the physical 
constants of the different materials are given in 
table 2. In choosing the constants, the speed of 
sound in infinite medium was used for c’ and the 
dynamic compressive vield strength was used for /” 
(In liquids there is only one speed of sound that can 
be used fore.) Dynamic rather than static strength 
values must be used for /£’ because the loading time 
is of the order of several microseconds. The strength 
that should be used for /” is the dynamic vield 
strength in shear. Unfortunately, the amount of 
work that has been done in determining the dynamic 
vield strengths of materials is very limited. Whiffin 
\7} has determined the dynamic vield strength in 
compression for a number of materials. In choosing 
the numerical constant for eq (21), the dynamie 
vield strength in compression was used for /” in- 
stead of the dynamic yield strength in shear. Some 
justification for this substitution may be found in 
the Von Mises strain energy theory according Lo 
which the tensile elastie limit, 
the vield strength in shear. 


.Is43 7, where +, is 


2.6. Test of the Equation 


For pit-depth calculations, eq (20) mav be 


in the form 


put 


where V, is the intercept velocity for the particular 
liquid-drop-solid-target combination being used and 
where V, is given by eq (21 Pit-depth-versus- 
velocity data obtained in another laboratory? were 
test the equations. In most the 


Cuses 


used to 


238 








TABLE 2. Physical Constants 
Sound speed, c Acoustic imped- Statice tensile | Dynamic compres- 
Metal Density, p infinite medium nce, vield sive yield 
strength, Y strength, - 
gicm cm/sec a/sec-cm pst dicm 
- ” £2,025 7. 239X105 
1100-0 alum im 2.713 6. 318X110 1.71410 nM) b 1. 18310 
2024-0 aluminum 2.768 1 ¢ ox10 1. 763 K 10 f 12. 625 «2? 35010 
Copy innealed electrolytic tough 
pit &. G2 14 691K10 4. 184 100 f a7 2.39410 
Lead, chen i 11. 3437 12. 27710 2. 583 X10 bh 4.324 10° 
Lead nt 2. 128 X10 2. 414x100 
St cold 1 i 7. 859 5. T8610 4. 54710 27410 
Iro st) e 5. 85010 1. SOR 10 H60*« 10 
Zire 7.14 4. 17010 2.977 X10 
Nlercury 13. 546 e1.451K10 1. O66 & 10° 
Water, 2 ( 0. 99707 1. 497 «10 0. 1493 10° 
Data from Aluminum Co. of America f Measured in NBS Engineering Mechanics Section by Lafayetts 
Data from Ilandbook of Chemistry and Physics kK. Irwin 
Data from Metals Handbook, 1948 « Data of A. C. Whittlin, communicated by letter. 
1 Nieasured in N BS Sound Section by Carroll Tschiegg Data of A. C. Whillin [7] 
e Data from L. Bergmann {12}. ; i Extrapolated from graph, figure 7e, given by Krafft and Sullivan [9 


physical properties of the metals used in the experi- 
ments were not determined. The physical con- 
stants used for these materials in eqs (21) and (22) 
are given in table 2. 


a. Test of the Equation With Mercury Drops and Pure Metal 
Targets 


Pit-depth-versus-velocity curves calculated by 
eq (22) for collisions of 0.10- and 0.15-em mercury 
drops against copper are shown in figure 8. Experi- 
mental points for the two drop sizes are indicated in 
this figure with circles and with triangles, respec- 
tively. The material used in these experiments? was 
deseribed pure copper. Physical constants for 
electrolytic copper were used in the calculations. 
The intercept velocity, Vj, was calculated by use 
of eq (21). It ean be seen that both the slope and 
intercept of the calculated curves are in good agree- 
ment with the experimental data. 
Pit-depth-versus-velocity curves calculated by 
2) for collisions of 0.10-, 0.15-, and 0.285- 
em mercury drops against 1100-0 aluminum are 
shown in figure 9. Experimental points for the 
three drop sizes that were used are indicated in this 
figure with cireles, triangles, and squares, respec- 
tively. The intercept velocity, V, was calculated 
by use of eq (21 The value of the dynamic com- 
pressive vield strength that was used for /” in 
(21) was calculated by means of a formula given 
by Whiffin for the Duralumins. The static 
vield strength of the 1100-0 aluminum used was not 
known.? For the purpose of the calculation of £’” it 
was taken to be 5,000 psi (table 2). The speed of 
sound in infinite medium, determined for a piece 
of 1100-0 aluminum of lower vield strength by 
Tschiegg (table 2) was used for ¢’ and in computing 
There is more scatter in the experimental data 
for the 0.15-em drop size than for the other two 
drop SIZES. Seatter could be caused by Variation 
of the drop size from the reported value; it could 
also be caused by use of target plates cut from sheet 
material of different vield strength. The data 
both the Q.15-em and 0.285-em drop sizes would fit 
the curves better if the intercept velocity were 
1.110" em/see rather than 1.26 10* em/sec, which 


ibs 


use of eq (2 


eq 


‘ 


for 


was the value found by use of eq (21). The experi- 
menters report [8] that the target plates were cut 
from 1100-aluminum sheets having thicknesses in 
the range of 0.185 to 0.210 in.; hence, they were not 
all cut from the same sheet stock and, consequently, 
some variation in yield strength, which would affect 
the intereept velocity directly, could have been 
present. In general, however, there is good agree- 
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FIGURE Y Calculated curves for collisions of mercury drops of 
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The pit-depth-versus-velocity curves calculated 
by use of eq (22) for collisions of 0.10-em and for 
0.285-cem mercury drops against lead are shown in 
figure 10. The experimental points for the two drop 
sizes are indicated with circles and with triangles, 
respectively. The intercept velocity, Vj, was cal- 
culated by use of eq (21). The value of the dynamic 
compressive vield strength used for /’ in eq (21) 
was that given for 99.97 percent lead by Whiffin 
(table 2). The value of the speed of sound in infinite 
medium for pig lead measured by Tschiegg (table 2) 
was used for ¢’ and in computing 2’. The calculated 
curves are a good fit to the data at the lower veloci- 
ties; at the higher velocities the depth of the pits is 
greater than that predicted by eq (22). The pits 
produced in lead by mercury drops are characterized 
by a heavier lip of metal around the mouth of the 
crater than is observed on pits in 1100-0 aluminum or 
on pits in copper. See figures 1, 2, and 3. It ap- 
pears that for a metal that is as soft as lead the flow 
of a liquid-drop projectile at high collision velocities 
drags a notable amount of metal from the bottom 
of the pit and piles it up at the mouth of the crater. 
This extra mode of pit formation was neglected in 
the development of eq and, therefore, this 
equation does not fully account for the depth of 
pits formed at very high velocities in metals as soft 
ibs lead. 


(22) 


b. Test of the Equation With Mercury Drops and Alloy Metal 
Targets 


Experimental pit-depth-versus-velocity data were 


obtained for collisions of 0.15-em mercury drops 
with steel targets in another laboratory? It was 
reported that the target material was originally 
4H cold-rolled steel. After annealing above the 


alpha temperature it developed a case hardening 
The case hardening was chipped off and the target 
plates were polished. It was reported that the 
resulting material had an average Rockwell E hard- 
ness of 90.96 and Brinell hardness of 90.98. which 
they stated corresponds to an ultimate tensile 
strength of 47,200 psi. 

In order to caleulate the intercept velocity by 
means of eq (21) it is necessary to know the dvnamiec 
compressive vield strength the target metal. 
The manganese and carbon contents of steels strongly 
affect their dynamic vield strength; the manganese 
content is particularly important if it is low, for 
example, in the range from zero to 0.50 percent. A 
chemical analysis of the manganese and = carbon 
content of one of the steel target plates was made 
in the NBS Analytical Chemistry Section. It 
vealed that the steel contained 0.70 percent man- 
ganese and 0.21 percent carbon. The value of the 
dynamic compressive vield strength used for /’ in 
eq (21) was found by extrapolating a curve given by 
Krafft and Sullivan [9] for a steel containing 0.98 
percent manganese and 0.22 percent carbon that was 
loaded in compression. The speed of sound in 
infinite medium determined for cold-rolled steel by 
Tschiegg (table 2) was used for ¢’ and in com- 
puting 2’ 


of 
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The pit-depth-versus-velocity curve calculated by 
use of eqs (21) and (22) for collisions of 0.15-cm 


mercury drops against steel targets is shown in 
figure 11 along with the experimental points. It 


would appear ‘that the intercept velocity is about 
correct. However, a straight line that would fit 
the experimental data would have a slope consider- 
ably than that of the calculated curve. It 
would appear from this evidence that eq (22) does 
not predict pit depth with exactness for pits formed 
in targets of alloy steel. 
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Che pit-depth-versus-velocity curves calculated 


by use of eqs (21) and (22) for collisions of 0.10-cm 
and 0.20-cm mercury drops against 2024—0 aluminum 
targets are shown in figure 12 along with experimental 
points obtained elsewhere? The experimental points 
for the two drop sizes are indicated with circles and 
triangles, respectively. The targets were prepared 
at the NBS. The average static yield strength 
(0.2% offset) of the metal used was found to be 12,625 
psi from measurements made by Irwin (table 
The dynamic compressive vield strength used for 
E’ in eq (21) for this material was calculated from a 
ratio of the static to the dynamic vield strength of 
2.7. This ratio was estimated by Whiffin [10] for the 
2024-0 aluminum that was used for the targets. 
The speed of sound in infinite medium determined 
for 2024—0 aluminum by Tschiegg (table 2) was used 
for ec’ and in computing The calculated curves 
differ from a best-fit line drawn through the experi- 
mental points in intercept but not in slope. If the 
intercept velocity had been 2x10 em/see rather than 
2.421x10* em see there would have been relatively 
good fit between the calculated curves and the 
experimental points. 

In summarizing, it can be said that eqs (21) and 
(22) have not given calculated curves that show as 
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A, calculated curve for 0.10-cm drop 
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good agreement with experimentally determined 
points for collisions of target plates of alloy metals 
with mereury drops as that obtained for collisions 
of plates of pure metals with mercury drops. 


c. Tests of the Equation With Waterdrops 


Test results with w aterdrops of constant size are 
available for some of the metals used in the experi- 
ments with mercury drops. 

The pit-depth-versus-velocity curve calculated by 
of (21) and (22) for collisions of 0.2-cm 
waterdrops against targets of annealed electrolytic 
tough pitch copper is shown in figure 13 along with 
experimental points obtained elsewhere. The tar- 
gets were prepared at the NBS. The static vield 


strength (0.2¢, offset) of the metal used was deter- 


use eqs 


mined by Irwin; the average value was 3,975 pst. 
The dynamic compressive vield strength used for 
FE’ in eq (21) was that given by Whiffin (table 2) 


copper. The speed of sound in in- 
determined for electrolytic tough 
pitch copper by Tschiegg (table 2) was used for e’ 
and in computing 2’. There is quite a bit of scatter 
in the data plotted in figure 13. However, both 
the slope and intercept of the calculated line are in 
relatively good agreement with these data. 

The pit-depth-versus-velocity curve calculated by 
of (21) and (22) for collisions of 0.2-em 


waterdrops against targets of 2024-0 aluminum 


for elect roly tic 
finite medium 


use eqs 


is 








shown in figure 14 with experimental points obtained 
elsewhere The targets were prepared at the NBS. 
They were made of the same 2024 aluminum as the 
targets that were used for collisions with mercury 
drops 2.6b and fig. 12), and the same 
values of /’, ¢’, and were used in calculating the 
pit-depth-versus-velocity curve. It can be 
from figure 14 that the calculated line fits the experi- 
mental points fairly well. 


see Sec. 
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The pit-depth-versus-velocity curve caleulate dl by 
21) and (22) for collisions of O.2-em 
waterdrops against leac targets is shown in figure 15 
The experimental points shown in the graph were 
obtained elsewhere The kind of lead that was 
used for the targets is not known. The dynamic 
vield strength used for /’’ was that given by Whiflin 
(table 2) for 09.97 percent pure lead. The speed of 
sound in infinite medium for pig lead determined by 
Tschiegg was used for e’ 

The amount of agreement has found 
between the ecaleulated and experimental 
points for collisions of both mercury drops and water- 
drops with targets of copper, 2024-0 aluminum, and 
lead miay 
usefulness 


Use of eqs 


and in computing 
that 
curves 


been 


be favorable evidence for the possible 
of eq (22) in predicting corresponding 
velocities-for-equal-pit-depth for collisions of metal 
target: plates with waterdrops and with mereury 


drops. 


2.7. Corresponding Velocities-for-Equal-Pit-Depth 


From the standpoint of the problem of high-speed 
rain erosion, the practical value of an equation giving 
the depth of damage pits as a function of lMpmgee- 
ment velocity is its use to extrapolate from = the 
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known depth of pit produced by collision of a solid 
target with a drop of high-density liquid at a rela- 
tively low impingement velocity to the depth of pit 
that would be produced on the same target material 
by collision with a drop of low-density liquid, such 
as water, at a very much higher impingement ve- 
locity. The purpose of carrving out this extrapola- 
tion is to bypass the necessity of firing target mate- 
rials at the extremely high velocities for which test 
results are desired for collisions of solids with water- 
drops. 

If drops of a liquid .1 having diameter d@, collide 
with a solid material having acoustic impedance 
and dynamic vield strength /’, eq (22) predicts that 
the pit depth will be 
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where Vy,4 is the intercept velocity calculated with 
use of eq (21) and with use of the density and acoustic 
impedance of liquid A for p and 2, respectively. The 
depth of pits resulting from collisions of drops of 
liquid B with the same solid material will be 


eo! 2d yz B , , 9 
OR 9 (\ Bo J al, (24) 
} 


Cea\<pte 


where V,, is the intercept velocity calculated with 
use of eq (21) and with use of the density and acoustic 
impedance of liquid B for p and 2, respectively. 

If the de ‘pth of pits sindened by collision of the 
solid target with drops of liquid A is to be the same 
as that produced by collision of the solid target with 
drops of liquid B, then, by equating 64’ and 6,’ given 


by eqs (23) and (24 
- dpZn€,(Z. 2! - 
Vig seematat ey, 
U4-ACRB\<-Bp TT ) 
d pz pes oA T See 2") 
2 a a PF an » 
dst iCp\<pB 7 ~ (Pppe < ies 
OR’ (+1 +’ 
19k a . # (25) 
Poe 1/9 _ 
Pa€ « ss 
All quantities are expressed in egs units. The dy- 


namic compressive vield strength must be used for /’ 
and the sound speed in infinite medium must be 
used for ¢’ and in computing in making calcula- 
tions with this equation. 


2.8. Applicability of the Damage Equation 


It was remarked in section 2.2 that the extent and 
tvpe of damage produced on solids as a result of col- 
lision with liquid drops depends strongly on the 
characteristic properties of the solid, and that if the 
damage is to be described by an equation, it will be 
necessary to develop separate equations for the 
damage produced on solid materials that have widely 
different: properties. It was pointed out that the 
simplest case to consider is that of materials such as 
the soft and medium hard metals that undergo 
permanent plastic flow without fracture as a result 
of collision with liquid drops. Equation (22) has 
been developed to deseribe the depth of pits that will 
be produced in materials of this kind as a result of 
such collisions. It cannot validly be applied beyond 
the limits of the model on which it was constructed. 

Kqui ation (22) is restricted to dams we associated 
with the impact pressure that results when a solid 
target collides with a liquid drop at high speed. 
The dimensionless Reynolds number, Weber num- 
ber, and surface tension quotient, which are associ- 
ated with the radial flow of the drop, were deleted 
from eq (7) to produce eq (8), and eq (22) contains 
only the dimensionless quotients that appear In eq 
(8S). In view of this restriction on eq (22) it can be 


expected that there will be limits on its application 
even to the soft and medium hard metals. 

In discussing the agreement of the calculated 
curve with the experimental points for lead target 
plates it was pointed out that for metals as soft as 
lead at very high collision velocities the radial flow 
of a mercury drop drags metal up the walls of the 
damage pit and piles it up at the mouth of the crater. 
This mode of forming, or of deepening, a pit was 
not considered in the development of eq (22), and 
eq (22) will not adequately describe pits that were 
produced either wholly or appreciably by this mech- 
anism. It is noteworthy that where this additional 
pit-forming mechanism does not operate (in the low- 
velocity range), eq (22) does adequately deseribe 
the pits that are formed in lead. 

Two other points should be mentioned with regard 
to the valid use of eq (22) for determining the depth 
of pits that will be produced in collisions between 
target plates of the medium hard metals (to which it 
applies) and liquid drops. One is the mode of 
mounting the plate of solid material with which 
the liquid drops collide; the other is the thickness 
of this plate. The reverse side of the plate must 
be maintained as a free surface for otherwise the 
core of solid material through the metal plate under 
the contact or collision area cannot move freely with 
respect to the remainder of the plate, which is an 
essential of the model on which eq (22) is based. 
The bulge on the reverse side of the 1100-aluminum 
plate shown in picture 6 of figure 2 appears to indi- 
cate that this condition was sufficiently realized in 
the mounting of the metal plates that were used for 
the mercury drop experiments.? If this condition is 
not realized, eq (22) cannot be applied to the pit 
depths that are obtained. The thickness of the plate 
should not be less than 1.4 to 2.0 times nor greater 
than 4 times the diameter of the impinging 
liquid drops. If the plate is too thin, the model on 
which eq (22) is based will break down because the 
yd will bend as a unit under the collision blow. 
If the plate is too thick, the spreading of the com- 
wie A Wave on passing through it may cease to 
be negligible. 


to 5 


Equation (22) was constructed on the assumption 
of permanent plastic flow of the material of the 
target plate. The dimensionless quotient ¢’//” of 
eq (7), which represents the resilience of the plate 
material, has been neglected. Consequently, eq (22) 
cannot validly be applied to determine the depth of 
pits that will form in highly resilient materials such 
as the polymers and rubbers when the collide with 
liquid drops. In the case of such materials the core 
of target material depressed as a result’ of the 
collision tends to spring back into its original posi- 
tion; the permanent damage mark that remains 
after the collision is more nearly a circular cut than 

pit. 

Equation (22) has been constructed using a very 
simple model. To determine how the Revnolds 
number and the Weber number should be introduced 
into eq (22), pit-depth-versus-velocity data should 
be obtained on a single metal for the condition that 
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the viscosity and surface tension of the crop liquid 
are gradually changed (such as by the use of glyeerol- 


water solutions). To determine how the sound 
speed, c, of the liquid should appear in the expression 
for the intercept velocity, data should be obtained 


using the same metals that have been used already 


and a liquid that has a sound speed different from 
those of water and mercury, which are nearly the 
same. Chloroform has a sound speed about 68 per- 
cent that of water. The sound speed of water can 


be increased by adding sodium chloride. The speed 
of sound in a 20 percent sodium chloride solution is 
about 13 percent greater than that in fresh water and 
it is indicated that increase in the salt concentration 
will further increase the sound speed. 

Finally, more — pit-depth-versus-velocitv data 
should be obtained for collision of both mercury 
drops and waterdrops against target metals whose 
dvnamic compressive yield strengths are known. 
Whiffin [7] has determined the dynamic compressive 
vield strength of standard silver (7.50% Cu, 92.5% 
Ag), electrolytic copper, 99.97) percent lead jand 
Armeo iron (0.0165) ©, 0.006°, Si, 0.0175) 3, 
0.0035, P. and 0.03007 Mn). In the process of 
verifving eq (22) it is important to use target ma- 
terials whose dynamic compressive yield strengths 
are known, and the same target material should be 
used throughout where drops of different liquids are 
used. If eq (22) is fullw verified, it can be used to 
determine the dynamic compressive vield strength 
of other metals 


3. Collisions Between Metal Plates and 
Flowing Metal Spheres 


VanF leet, and Whited [11] fired spheres 
le ad. aluminum, and iron against 


Partridge, 


of zine, tin, copper, 

targets of the same material at = collision veloc- 
ities up to 2410" em/see (7,900 ft/sec) and main- 
tained conditions such that the spherical pellet 


lost no mass before striking the target. They found 
that the penetration varied linearly with velocity for 
the materials used and the velocity range investi- 
gated. They reported that the pellets flowed during 
the collisions 
In the light of the fact that the pellets flowed 
during the collisions, it is reasonable to suppose that 
the pit-cdepth-versus-velocity equation that was 
developed for collisions of metal targets with liquid 
drops, eqs (21) and (22), should apply to this case 
without modification. When the material of the 
target is the same as the material of the projectile, 
eqs (21) and (22) simplify to 
= 3.6 d (V—V,) /c 26 
and 
V,.=38 E’/2 (27 


Some of the pit-depth-versus-velocity data [11] are 
plotted in figures 16, 17, 18, 19, and 20 along with 


curves calculated by use of eqs (26) and (27). The 
physical properties of the metals used by Partridge 
et al. [11] were not determined. The physical 
constants used for these materials in eqs (26) and 
(27) are given in table 2. 


3.1. Collisions of Iron Spheres Against Iron Targets 


It is assumed that soft, relatively pure iron was 
used in these experiments [11] because the spheri- 
cal pellets fired were made by placing fragments 
of the metal in the hemispherical cavities of a 
case-hardened steel tool and pressing the two 
sections together. The diameter of the pellets for 
which the experimental pit-depth-versus-velocit) 
data are given in figure 16 was 0.483 cm. The 
calculated curve is shown in figure 16 along with 
the experimental points which are indicated with 
circles. The observed pit depths for iron spheres 
colliding with iron targets are in good agreement 
with the calculated curve up to a collision velocity 
of 2410' em/see (7,900 ft/sec). This velocity is, 
however, less than half the speed of sound in infinite 
medium for iron (58.5 10° em/sec). 
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Figure 16. Calculated curve for collisions of O.483-cem tron 


sphe res against tron targets 


observed depth for 0.483-cem sphere 


3.2. Collisions of Aluminum Spheres With Aluminum 
Targets 


Values 8 — ‘al constants of 1100-0 aluminum 
were used 1 (26) and (27). The value of the 
dynamic Posteo sr vield strength used for 
in eq (27) was computed from a ratio of the dynamic 


egs 


to the static vield strength of 4: this is a rough 
ratio estimated by Whiffin [10] for an 1100-0 
aluminum having an average static vield strength 
of 2,625 pst. 


Very few experimental pit depths were given by 
Partridge et al. [11] for collisions of 0.483-cem alumi- 
num spheres with aluminum targets. The data 
available contain quite a bit of seatter. They are 
plotted in figure 17 along with the curve calculated 
using eqs (26) and (27). There is reasonably good 
agreement with the calculated curve. The highest 
collision velocity for which a pit depth was re ported 
is approximately one-third of the speed of sound in 


aluminum. 
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Figure 17. Calculated curve for collisions of 0.483-cem 
aluminum spheres against aluminum targets. 


observed depth for 0.483-cm sphere 


3.3. Collisions of Lead Spheres Against Lead Targets 


Whiffin [7] determined the dynamic compressive 
vield strength of lead that was 99.97 per cent pure 
and this value of the dynamic compressive vield 
strength was used for £’ in eq (27). The sound 
speed in infinite medium was determined by Tschiegg 
(table 2) both for chemical lead and for pig lead. 
The theoretical curve computed using the sound 
speed for pig lead is shown in figure 18 along with 
the pit-depth-versus-velocity data obtained by 
Partridge et al. [15] for 0.483-cm lead spheres 
impinging against lead targets. 
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Figure 18. Calculated curve for collisions of *\9-in. lead 


sphe res against le ad large ts. 


«-ln sphere 


The experimental pit depths found at the high 
velocities used are all greater than eq (26) would 
predict. No data were obtained in the low velocity 
range where the calculated curve was found to fit 
the data for collisions of lead target plates with 
mereury drops. See section 2.6a, where the pos- 
sibility that two mechanisms contribute to the 
formation of pits in lead at high collision velocities 
is discussed. 


| 
| 
| 


The speed of sound in infinite medium for pig lead 
is marked with a vertical dotted line in figure 18. It 
would appear from the empirical data at the highest 
velocities used that the experimental pit-depth- 
versus-velocity curve for lead is flattening at veloci- 
ties above the speed of sound in lead. 

[t appears from the data that have been presented 
in the preceding sections that, at collision velocities 
for which an originally solid target remains solid dur- 
ing the collision, pit depth is a straight-line function 
of collision velocity if the rear face of the target plate 
is a free surface. This appears to be the case regard- 
less of whether the projectile was originally solid and 
remains solid during the collision (low-speed solid-to- 
solid collisions), whether the projectile was originally 
solid but flows during or as a result of the collision 
(high-speed solid-to-solid collisions), or whether the 
projectile was originally a liquid drop (liquid-to-solid 
collisions). At extremely high (meteoric) collision 
velocities it appears that an originally solid target 
may be expected to behave like a liquid during the 
collision because of the enormous impact pressure 
developed. With regard to collisions of this kind 
Opik [13] has stated that the aerodynamic (1/2 pV*) 
pressure at the penetration of a meteor into rock ts 
more than 1,000 times the plastic limit of steel; he 
has hypothesized that all solid materials under such 
pressures must behave like liquids and that the 
problem of meteor impact is the case of the impact 
of a liquid drop of one density into a liquid medium 
of a different density. Under such conditions the 
penetration or pit depth may become independent 
of the collision velocity, or may become essentially 
so, and the pit-depth-versus-velocity curve may ap- 
proach a horizontal line parallel to the velocity axis. 
A trend in this direction may be indicated by the 
flattening of the experimental pit-depth-versus- 
velocity curve for lead. Whenever this condition is 
realized, either partially or completely, the pit- 
depth-versus-velocity equation developed in_ this 
paper will no longer be applicable. 


3.4. Collision of Copper Spheres Against Copper 
Targets 


The data for impingement of 0.483-cm copper 
spheres against copper targets [11] are shown in 
figure 19 along with the theoretical curve calculated 
by eqs (26) and (27) using physical constants for 
electrolytic copper (table 2). 

The empirical pit-depth-versus-velocity curve for 
copper appears as though it may be flattening at the 
highest collision velocities for which data were ob- 
tained, as was found in the case of lead. If this ts 
indeed true, it raises the question as to what deter- 
mines the flattening of the curve. In the case of 
lead it happened at collision velocities equal to the 
speed of sound in lead. If it is happening for copper 
at velocities of 16 10* to 18 10* em/see, it ts 
happening at velocities that are about half the speed 
of sound in copper. The speed of sound in infinite 
medium for electrolytic tough pitch copper is 46.91 

10' cm/sec. It does not, furthermore, seem to be 
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related to the melting point of the target material 
because copper has a high melting point and lead 
has a very low one. A possible explanation is the 
relative susceptibility of the target material to flow 
by translational slip of the target atoms on their 
lattice planes at high rates of loading. With regard 
to the use of copper as a liner for shaped charges, 
Rinehart and Pearson [14] state that copper flows 
readily at high rates of loading. 
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3.5. Collisions of Zinc Spheres Against Zinc Targets 


No value of the dynamic compressive  vield 
strength of zine was obtainable. One of the possible 
uses of the pit-depth-versus-velocity 
ported in this paper is the calculation of the dynamic 
compressive vield) strength from pit-depth-versus- 
velocity data. By trial it was found that to obtain 
an acceptable value of the intercept velocity for zine 
with the pit-depth-versus-velocity data of Partridge 
et al 11], the value of the dvnamic COMpressive 
vield strength for zine would be 1.546 10° dem. 
The experimental pit-depth-versus-velocity data for 
zine and the line calculated from eqs (26) and 
using this value of the dynamic compressive vield 
strength and the physical constants given in table 
2 are shown in figure 20.) The agreement between 
the experimental points and the theoretical curve is 
reasonably The two points obtained at the 
highest velocities used seem to show a flattening of 
the curve. More experimental data are needed to 
determine whether this ts the case 
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3.6. Collisions of Tin Spheres Against Tin Targets 
al. 


data 


[11] also obtained pit-depth- 
versus-velocit, for collisions of tin spheres 
against tin targets. The tin appeared to crack and 
break off rather than to flow plastically around the 
crater. This interesting fact recalls the similarity 
in behavior of materials at low temperature and at 
high rates of loading. At low temperature tin tends 
to exist in the brittle gray tin form. It is doubtful 
whether the pit-depth-versus-velocity equation can 
be validly applied to tin because of the e1 erge\ that is 


Partridge et 


diverted from pit formation into crack formation 
see section 2.S 
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Earth Currents Near a Monopole Antenna With 
Symmetrical Top Loading 
James R. Wait 


I}xpressions for the fields are developed for a vertical ground-based monopole with a cone or disk located 
at the top of the antenna to simulate umbrella top loading. The current distribution on the structure is 
assumed. Using spherical-wave functions, the magnetic-field distribution on the ground plane near the base 
of the antenna is computed and illustrated by graphs. For the case where the antenna is electrically small, 
the currents flowing on the cone or disk are shown to contribute onty slightly to the total field 


1. Introduction 


Transmitting antennas at low radiofrequencies usually take the form of a vertical current- 
carrving wire loaded by extensive overhead wires which sometimes take the shape of an um- 
brella. At the risk of oversimplification, a model for field calculations might be a vertical 
monopole with the umbrella being represented by a continuous cone. In other words, any 
asvmmetry about the vertical axis is neglected. It is the purpose of the present note to inves- 
tigate the expected influence of the umbrella or cone loading on the fields produced by the 


antenna with an assumed current distribution. 


2. List of Principal Symbols 


. 4, and @,—spherical coordinates of a variable point on the conical sheet. 
a—distance from the apex to the rim of the cone. 
J,—current at the base of the vertical member. 


F(=)=current distribution along the vertical member. 
h=height of the antenna or length of the vertical member. 
II! 2 component of the Hertz vector of the currents on the vertical member. 


K=the slope of the current distribution on the vertical member (assumed 
constant 
f(r) =the current distribution along a generator of the cone. 
A,= coefficient for the power series expansion of /(7)). 
a variable point on the vertical member (the integration variable). 
P, QD, = evlindrical cocrdinates of the observation point af 
/7;—tangential magnetic field on the ground plane resulting from the currents 
on the vertical member. 
IS and WS—the p and 2 components of the Hertz vector of the currents on the cone. 
R—distance from the variable point ( on the cone to the observation 
point 2. 
r, 0, 2=spherical coordinates of the point 7’. 
Py (cos #)=associated Legendre polynomial of argument cos @ and orders n and m. 
yntkr) =spherical Bessel function of the first type of argument ér and argument n. 
h\? (kr) =spherical Hankel function of the second kind of argument kr and order vn. 
T(r, 7) =yn (kr )h& (kr). 
/7¢—tangential magnetic field on the ground plane resulting from the currents 
on the disk or the cone of angle @,==90 
Ax(ro), Balto), Cr(%o) = coefficients for the spherical-harmonic expansion of 77/7. 
ry and @)=spherical coordinates of a variable point on the ground plane for the 
origin at O in figure 1, ry=(p?+-A?) and @—aretan (p/h). 
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j' QHi/I, . . earth current (dimensionless) on the ground plane resulting from the 
currents on the vertical member. 
j all? 1 earth current (dimensionless) on the ground plane resulting from the 


currents on the disk. 


7, @—new spherical coordinates for describing the radiation pattern (se 
fig. 3) 
AS), AOS), and JCS) are evlindrical Bessel functions of the first type of orders 1, 2, and 3 of 
the argument S 


3. Some Preliminary Considerations 


The top-loaded antenna in its idealized form is shown in figure 1. It consists of a vertical 
current-carrying wire extending from the point QO, at 20, along the axis of a evlindriecal 
coordinate system to the ground surface at z=A. The point O is also the apex of a cone 


which is generated by a line of length, a, at an angle @, to the 2 axis. With respect to a spherical 
coordinate system with the origin at O, the coordinates of the bottom rim of the cone are (a, 
4,. @) and those of a variable pot (Jon the cone are (7;, G,, D,). The upward current on the 
central wire filament is taken to be /,/(2) where /, is the current at the apex, Le., F(0)=1, 
and /,/(h) is the current at the base. The outward current on an element dr,d¢, on the 
conical surface is then (/)/2)d¢@f(7,) where /(0)=1 and f(a)=0. Because of symmetry the 


current on the conical sheet has only a component in the direction of increasing 74. 
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The current distributions, characterized by the functions F(z) and f(r), in general, have 
a very complicated dependence on the antenna dimensions. In the absence of a top loading 
and where the vertical member can be approximated by a circular cylinder, current distribu- 
tions can be expressed in terms of an expansion parameter which is proportional to log (h/27c) 
where ¢ is the radius of the conductor [1].!. When A is less than one-half wavelength, the 
distribution is very close to being sinusoidal so long as / is greater than about 100 ¢. In fact, 
at very low frequencies where h becomes small compared with the wavelength, the current is 
then almost a linear function of height. It would not be expected that the cone would modify 


this conclusion. Therefore 


F(z) =I,1—K 2/h), (] 


where A is a constant which depends on the extent of the top loading. For example, in the 
unloaded case, A=1, and in the ideally loaded case, A= QO, corresponding to a constant current 
on the vertical member. In most applications where top loading is employed it is known 
on experimental grounds [3] that AO. This is particularly so for high-power installations 
at low radiofrequencies where the effective dimension of @ is comparable to or greater than 
h, and @, is near 90°. 

Information regarding the current distribution function f(7,) for the conical sheet is 
apparently not available in the literature. However, if each individual wire making up the 
cone acts independently, it could be expected that the total current flowing towards the rim 
of the cone would vary in a linear manner with 7, if a is small compared to the wavelength 
On the other hand, when the cone degenerates into a disk (1.e., 6,~90°) and, if A is somewhat 
smaller than a, it can be expected the total current flowing towards the rim will vary with 


the square of 7;. This conclusion is based on the idea that the shunt displacement current 
densitv between the circular disk and the ground is a constant if fringing effects at the edges 
are neglected. A plausible assumed distribution would be 
r q 

f(7;) ( ) Ay, (‘)) 

iid | a 
where the summation could be over g=0 and 1 and Ayp=1, A, | for the linear current 
approximation; whereas g=0 and 2 with Agp=1 and A, | for the squared current approxima- 
tion. It would be possible to employ more elaborate forms of f(7;) by utilizing more than two 
values of g. A prime requirement is always 

> 3) 


since f a (). 


4. Field of the Vertical Member 


> 


The Hertz vector II’ at the point P(p, 0, 2) due to the current flowing only on the vertical 
member is given by II*=(0, 0, IID), where 


1 f y =P ALC ; a : et] F(z’)dz’+ i eS te rer F(z’)dz’ }; 
for 1we, J0 [(z . 4 p*| ‘ (2h E yy Te. 


when /£&=2m7 wavelength and @=8.85410~-". The fields are then given by 
; O*IL: ; O° , oll, 
ID ©» Et=(P+—,)l:, and Hj=—iew =~. 5 
" Opoz 92 ; ? Op 
Figures in brackets indicate the literature references at the end of this paper. 
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The second integral in eq (4) can be identified with the image of the vertical member below 
It can be verified that “> =0 at h, which is the required boundary condi- 


the plane paz. 
tion for the ideally conducting ground. 
The tangential magnetic field, 77,, on the ground plane is then explicitly given by 
/, O ™ Fis’ exp ikl | 2’ —h)?+4 p?}’? P : 
Hi=— , Th dz’, 6 
2x Op J [(2 h)?+- p?|’* 


Considering the important case where Bh- l and /(2’)=1, it follows that 


IT; _— a [ 4 /\9 . a 9 uk 9 dz’ 
2a J [p” at. 2] Ps I” 
eke] h ‘ fh 
| 52 A) tikp tan & )} 


2 rp 


~ 


This can be expanded in a power series in (h/p), so that 
| 


; kp ) T terms containing ( m ). é ‘, ote. | S 
) p 


( nd | h “7, " h l l 
Hi~— ( ) (+ ike) (7) (5+: or ; 
The 


The first term of this expansion corresponds to the field of a hertzian dipole radiator. 
additional terms can be regarded corrections to account for the finite value of the antenna 
height and can be called multipole terms. Sufficiently close to the base, such that /p< <1, 
the following expression will suffice 

/ h 

% , qg 
5 : 


-7T p p” i? 


HH, =- 


which is a quasi-static approximation \ check on this result is obtained by noting that 


-7p 


lim 77; wat 


which could be predicted from Ampere’s law. 


5. Field of the Conical Currents 


From an argument based on SVinmetry, 1 can be seen that the Hertz vector If at P., due to 
The p and components, 


the currents flowing in the cone, has a vanishing @ component. 
These can be expressed as an integration over the 


however, are, in general, nonvanishing 


area of the cone. such that. 


[,/2m) sin @, (°"" “ e— ikk 
I] do dr, cos @ f(r 10 
: hor ie Ww e e R 
and 
1/22) cos 6, (2 a ec . 
I | lo dr, —,— fir), 11 
hor ie ‘. . . Ii 


ae 2rr, cos Qand cos Q= cos 6 cos 4 sin #@ sin @ cos o,. Now the factor 1 Risa 


where 2? =r 
solution of Laplace’s equation and can be written in terms of spherical harmonics, as is well 


known [2], 
. ‘ n—m)! ,,, 9 | ee a 
S S Pp Cos @ p= (cos @ COS I Oy 12) 
aed | a tele ts if 


PR —— tat * ! 
‘ ( i=0 i Mi). 


Where @ = 1, €,=2 1f (m0), and where the P? are associated Legendre polynomials. The factor 
R, on the other hand, is a solution of the wave equation and an expansion for it in terms 
Such an expansion does not seem to be available; 


f 


of spherical wave functions can be expected. 
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it can be obtained, somewhat heuristically, by noting that it should reduce to eq (12) for kr 
and #&r,; small, and should give rise to outgoing wave at infinity. Such an expansion is 


ikR n ! " ay (2) ay)? yp , 

é ; (n—my)! pnlkry hy (kr); r>r; 
- —ik >\ >) em » (2n-4 1) P” (cos 6) P™ (cos 6;) cos md, ; = (13) 
h n=0 m=0 (m-+mM)! pukry hy? (kr;); Ki Fr, 


where jn(kr;) is a spherical Bessel function and A,’ (kr) is a spherical Hankel function of the 
second kind. These are defined by 
, d "/sin ws 
gn(r) = (—1)"2r" ( ) ( ); (14) 
. rdr I 


and 


iz 


hj? (x) =171(—1)” ( — ) ( ; : ) (15) 


The integration, with respect to ¢,, in eq (10) and (11) requires the evaluation of 


| °or ikR 
TI, (r,0) ‘a ; R cos o1¢, (16) 
oT Jo i 
and 
l 20 | ikR 
1, (7,0) =i PR dd,, (17) 
= JO ‘ 


which can now be rewritten 


(2n+1) T,(r,,"3;7r>r 
1, (7,0) ik >> P} (cos 0)P4 (cos 6,) 4 5") ' : (18) 
fai n(n+1) T.(?, %)3 %) >? 
and 
T,.(%1, r);7>r 
1,70) ik S* 2n+1)P2 (cos 0)P2 (cos 0,) ath edith (19) 
” po | rf, Py) ly r 
where 7,,(7,. 7) =jn(kr he (kr). 
The circumferential field /7{ due to the conical currents is given by 
— oll; Olt: 
HS =ieqw ?— a (20) 
Oz Op 
: oll; Ol, . Oll: . oll: 
Hi Je w ? cos 6——— sin 6— -— sin 0— — cos 6 }5 (21) 
or roe or roe 
where 
: J,sin@, (* , 
Leow II | fir lL r@odr,, (29) 
4x Jo 
and 
I, cos 6, (** 
‘egwllS= : | fir) Ter Odr,. 23) 
4er Jo 
To satisfy the boundary conditions at the ground plane, z =A, the fields of the mirror image of 
the cone must be added to //‘. For these fields /,(7,@) must. be replaced by L(7’,8’) and 
I,(r,0) by + 12(r’,6’) where 0’ =M—8@, r’=[(2h—2)?+-p?]?, and cos M=[(r)?+ (r’)?—4h2|/2rr’. 


The preceding equations, along with eq (7) and (8), constitute the complete formal solution 
for the fields of a cone-loaded monopole fora specified current distribution on the strueture. 


6. Disk-Loaded Monopole 


When the angle @, approaches 90°, the cone degenerates into a disk. It is this case that 
will be considered in some detail here. Of particular interest is the behavior of the currents 
excited on the ground plane in the vicinity of the disk. In faet, for vlf antennas these near- 
zone currents have, associated with them, considerable ohmic losses due to the imperfect condue- 
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tivity of the soil [4]. The surface current density is by definition numerically equal to the 
tangential magnetic field /77 on the ground plane. For the case kr<C <1, it follows that 


, Lo . : 
Ht=~" S> [A, (7) Pha. (cos 9) + Br(to)Ph-1 (Cos 6) |Pi.(0) for ro<a (24) 
2r n l 
and 
HT ay Uno Pp : (COs 6.)P4(0) for r, >a, (25) 
27 Aa) 
where 
a f(r 7" 
A, (Yo —| vt dr,, (26) 
f (n+1)rq*? 
“@ f(r)rh 
BLO oS dy (27) 
: nr 
"e $fr,)r 
Cc, (Ta -- Bd + dr, (28) 
J 0 \li- L7G P , 
and 7)= (p?+-A?)'? and @=are tan (p/h 
Setting 
ae >A, (r; a)*, (29) 
it is a simple matter to arrive at 
l . (r)/a)*A 
A, (rm) =-— > 30) 
ry “G (n+1)(n+q7+!1 


— 
“e 
~ 
M 
— 
~ 
 —— 
=~ ~ 
~_— | 
— 
— 
=~ ~ 
=~ | 
— 
Senne 
— 


( ; a. ~, l ee | 9 

“7 To (n+q+1) (n+! 
Setting o1 l, A l and the other J1,’s equal to zero, the squared current distribution, 
Aw | ri/a*), is attained. Numerical values of the normalized current, } ally) 1, 
and = 7' alli/1, are shown plotted in figure 2 as a function of p/a for the squared current 


distribution on the disk and a constant current on the vertical member. The values of the 
ratio A/a are shown on the curves. The resultant current 7 on the ground plane which. is 
radially direeted is given simply by /*+,/° in the near-zone region where p is small compared 
to the wavelength. 

It is of interest to examine the radiation field such that r>>wavelength. The spherical 
Hankel functions can then be well approximated by the first term of their asymptotic expan- 
sion. The radiation field /7$ of the disk by itself is then given by 


ihe lye : an 2n-+ 1) =m - "a , ' om 
TT? — cos 6 >) : P' (cos 6) P' (0) qn kr) fry dry. 30 
Tl Tm LAeL | 


The radiation of the disk and its image is obtained by adding He to the preceding equation 
which is of the same form as —//f with r—@ replacing @, and r—2hk cos 6 replacing r. The 
summation in the preceding can be converted to a Bessel function by first starting with the 
addition theorem; 


’ ' 
. ° i WL). -_ " ° . 
exp ihr, sin 8 cos @ a € " (2n+-1)P™ (cos Q) f" 0) cos MO pn kr), (34) 
m=0 n=' Ji 770): 


252 


























Ss 
rm 


Normalized Earth Current jor j° 











0 ip 10 


Distance From Base _ P 


Radius of Disk fe) 





I: LGURI 2 ka th current distribution he lowad sk-loade a MOnO pote 


253 








and then multiplying each side by cos @ and integrating with respect to ¢@ from 0 to 27. This 


leads to the interesting relation 


: , (2n+1) a, 
iJ, (kr, sin #)=>3 Ph (cos €) P30) jntkr 35 
wae i 
There fore 
he Le ; ™ 
H] cos 6 fir JJ (kr; sin #dr, 350 
tr) J e,a8 
subject, of course, to kr >> 1. The squared-current assumption, f(7;)=1—(73/a*) is now intro 
duced in the above which fortunately enables the integration to be carried out to \ ield 
I kae~** - 
i cos @7(S), oF 
Lr? 
where 7S |] J(S J, SOS with S=fka sin 6. The field of disk plus its image can then 
be written 
1, _ Iokae =. = ~ 
H° +H >= cos 67 sin (kh cos 6) 7 (ka sin.@), Os 
— T/ 


in terms of the spherical coordinates 7, @ with the origin at the base of the vertical member as 
indicated in figure 3. For the usual case where ka sin @- |, it follows from the power series 
expansions of the Bessel functions -/J, and J, that 

ka sin 4 


Tika sin 6 g ) 34 


“=I 











—, 20 — - ~ 
_ 
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| 
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where terms in (ka sin @)*, ete. have been neglected. Noting that the radiation field for the 
vertical member in an 7, @ coordinate system can be written 


ik The 87 


2rr 


H: 


sin @, 10 
it follows that the resultant radiation field is given simply by 


miewontg (1S ee ' 


S 
neglecting terms of order (ka)?*. 


7. Conclusion 


The results indicate that the earth currents near a monopole antenna are not significantly 
modified by the presence of a symmetrical top loading. Therefore, in any calculation [4] of 
the energy absorbed in the ground system, it is usually sufficient to consider only the contribu- 
tion from the currents in the vertical members (downleads). Furthermore, the radiation pattern 
of a loaded monopole at low radiofrequencies is essentially a dipole pattern being modified to 
small order by the currents in the top-loading structure. The effect of unsymmetrical top 
loading is considered in a sequel to this paper [5 


The author thanks Mrs. A. Murphy for carrying out the calculations and W. W. Brown, 
A. D. Watt, and H. A. Wheeler for informative discussions. 
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Infrared Absorption Spectra of Some Cyclic 
Acetais of Sugars 


R. Stuart Tipson, Horace S. Isbell, and James E. Stewart’ 


Phe 


infrared absorption spectra of twenty-eight 
dene acetals of sugars have been recorded and analyzed 


l-methoxve thy lidene and ISOpPropy li 


The formulas of the compounds 


were grouped according to structure, and the absorption bands for each group were con- 


sidered in relation to the bands found for the other structurally related compounds 
bands were found for the unequivocal detection of the 
a ketopentofuranose, 
an aldohexopyranose or a ketohe xopyvranose, a& Uronic acid salt or lactone, 
Regardless of the fine structure of the dioxolane com 


tively to an aldopentofuranose o1 


furanose, 
keto-acid salt 


aleohol, or aw SUNT 


pounds, certain strong bands were observed, but the bands are 


No 
1,3-dioxolane ring attached respec- 
awn aldopentopy ranose, an wldohexo- 
a SULT 


not easily recognized from 


th Inspection of one or two spectra and are not characteristic of these compounds only 


1. Scope and Purpose of the Project 


For the purpose of discovering correlations that 
might be of value in structural analysis, this project 
was primarily undertaken with the object of provid- 
ing infrared spectrograms of sugar derivatives con- 
taining the 1,3-dioxolane ring (1). 


It seemed probable that the vibrations of the di- 
oxolane ring of such compounds would be relatively 
insensitive to changes In the rest. of the molecule: 
and that it should, accordingly, be possible to identify 
certain bands in their spectra as being characteristic 
of this ring Then spectra have, therefore, been re- 
corded; and the bands have been compiled and 
studied by statistical and comparative methods 


2. Compounds Investigated 


In the application of infrared absorption spectra 
to the determination of structure, it is necessary to 
have reference spectra for substances of known con- 


Presented, in part, before the Division of Carbohydrate Chemistry at the 
134th Meeting of t Amer n Chemical Society, at Chicago, Ill., on Septen 
her 9, 1958 Part of this work was supported by the Chemistry Branch of the 
Mtlice of Naval Rese 

* Present ire Beckn Instruments, It Fullerton, Cali 


figuration and conformation. The pyranose ring in 
simple carbohydrate derivatives is flexible; hence, 
for some pyranose derivatives, there is much uncer 
tainty concerning the conformation adopted by the 
ring. In contrast, the conformation of molecules 
containing fused or bridged rings * is more rigid; and 
consequently, less uncertaints exists as to the rela- 
live steric positions of the atoms and groups within 
their molecules. For this reason, attention has now 
been directed to the study of a group of evclic acetals 
(of evelie having bievelie and trievelte 
structures, 

Table 1 gives a list of the compounds, their code 
numbers [1] *° and properties, and an index to the 
spectrograms (the serial number of a compound is 
the same as the number of its spectrogram The 
spectra were measured in the 2- to 15-u region (5,000 
to 667 em7!; sodium chloride optics) and in the 15- 
$0-u region (667 250 em cesium bromide 
optics Inasmuch detailed study of all the 
bands is not possible at present, the spectrograms 
are given together with a discussion of (a) the strue- 
ture of the compounds and (b) some of the outstanding 


sugars ) 


to to 


as oa 


features of their spectra 
have two adjacent atoms in common, whereas bridged rin 


Fus 


onadjacent atoms in common 
Figures in brackets indicate the literature references at the end of tt paper 
Phe references for table 1 are given at the end of the table 
Code numbers are assigned according to a previously published fication 
t l!, for use with punched-card 1 hreiqguae 
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As a common structural feature, these compounds 
have the substituted 1 .3-dioxolane ring: 


where R—OCH, in the 1-methoxvethylidene acetals 
methyl orthoacetates), and R=CH, in the isopro- 
pylidene acetals. All of the compounds have at 
least one dioxolane ring; and, except for 1,2:5,6-di-O- 
isopropvlidene-p-mannitol (compound 28), this ring 
is fused to either a five- or a six-membered ring. 

The precise shape of the dioxolane ring is not 
known. The bond lengths and bond angles of the 
atoms comprising the ring are, however, compatible 
with a planar conformation. The abundance of 
compounds of this kind having structures explicable 
by a planar or nearly planar dioxolane ring would 
tend to support the hypothesis that the ring ordinarily 
has a planar conformation, 


3. Structures of the Compounds 
3.1. Methyl Orthoacetates 


The I-methoxvethvlidene acetals or 1,2-(methyl 
orthoacetate(s)), compounds | to 6, were prepared 
in investigations on the mechanisms 
of replacement reactions conducted over the course 
of many vears [2, 3, 4 They are acetals containing 
the 


the course of 


() OCH 


group fused to Cl and C2 of a pyranose ring. In 
both chair conformations,® the eis pair of hydroxyl 


sion ot the ontormatiors ot pvranoid ring sgiven ina previous 


14. F. Shafizadeh and M. L. Wolfrom, J. Am. Chem 77, 2568 (1955 
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groups at Cl and C2 of the parent sugars is equa- 
torial-axial. Spanning them witha planar dioxolane 
ring would introduce strain. Hence the pyranose 
ring or the dioxolane ring, or both, would be some- 
what distorted. The following depicts the skeletal 
structure of a possible sofa conformation: 


oe ~ 4 / 


For compounds 1, 4, 5, and 6, R 
pounds 2 and 3, R=H, and the structure is 
mirror image of that depicted. 

Although the pyranose ring normally assumes a 
chair conformation, spanning of Cl and C2 of a 
boat form of the pyranose ring with a planar dioxolane 
ring tends to place Cl, C2, C3, and the oxygen 
atom (of the pyranose ring) in a plane at an angle 
of 120° to the dioxolane ring. Two conformations 
of this type seem possible, depending on the positions 
of C4 and C5 of the boat. In one conformation, 
hereafter designated endo, they are directed toward 
the dioxolane ring, with C6 in an axial position; 
in the other conformation, designated these 
carbon atoms (C4 and C5) are directed outward 
from the dioxolane ring, with C6 in an equatorial 
position,’ The two skeletal structures are depicted 
in the following formulas: 


OAc: for com- 
the 


€L0, 





bo} 


In forming the 1,2-(methyl orthoacetate) from 
the poly-O-acetvlglvcosyl halide by the opposite-face 
mechanism |{6], the endo structure would arise by 
turning the nascent dioxolane ring (in the transition 


An end 
ubstituent es outs 





the V of a bieveliec structure, whereas 


’ 


ubstituent ties In 
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structure 
The ero 


state) inward toward C6; and the exo 
(Il) would arise by turning it outward. 
form would seem to be the more probable for the 
methyl orthoacetates, because the overall structure 
is more nearly planar and there is less crowding of 
the atoms 

All of the orthoacetates studied have essentially 
the same ring structures. A summary of the posi- 
tions of the reference groups in the ero-boat form 
Il) of each of these methyl 
vivel in table 2. The ester carbon atom 
with an asterisk) of the orthoacetates is asi mmetric: 


orthoacetates is 
marked 


there is no experimental basis for assigning configura- 
tion to it, but since the methoxyl group ts larger 
than the methyl group, it would tend to be directed 
away from the V of the bievelic structure. 

Decision between the possible conformations mndi- 
cated must await the outcome of further experimental 
work 


derivatives analogous to compounds 1 to 


[sopropy lide The 


7. 1,2-0-lsopropylidene-p-galactopyvranose 

Compounds 8 to 10 are isopropylidene acetals 
which, unlike the orthoacetates, have two dioxolane 
rings fused, on opposite sides, to a pyranose ring. 
The presence of two nearly planar dioxolane rings in 
this position is not compatible with the ero-boat 
conformation of the pyranose ring most probably 
present in the orthoacetates (compounds 1 to 6 
The two pairs of cis hydroxyl groups (at Cl and 
C2, and at C3 and C4, for compounds 8 and 9) of the 
parent each axial——equatorial in both 
chair-conformations;® and spanning them with a 
planar dioxolane ring would introduce strains (see 
IVa). A sofa conformation, in which the pyranose 
ring or the dioxolane rings (or both) may be some- 
what distorted, seems possible. One structure with 
very little strain is shown in formula IVb 


sugars are 


6 (having a methyl group in place of the methoxyl group PS . 

seem feasible but are not known at present The 1,2-0- 2 L > J 
isopropylidene-bp-galactopyranose considered in the next see- yr ¢ HAS HA 
tion is closely related; however, it differs in having the -R 

opposite configuration at Cl and C2. The 1,2:5,4-di-O- Yt ¢ , 

isopropylidene derivatives discussed in the next seetion are 


ilso related to the orthoacetates studied, but thev mav have 


t more rigid strueture because of the stabilizing effect of the 


SCOT d dioxolane ring 


3.2. O-Isopropylidene —Pyranose Acetals 


Compound 7 of table Ll has one dioxolane ring 
fused to a pvranose ring The latter could have a 
boat (probably Lo conformation: however, il might 
have a distorted sofa conformation in which, at 
C2, C3, and C4, both eroevelic bonds are approxi- 
mately quasi instead of axial equatorial. 
The conformation IT] that is equatorially 
y-substituted is depicted 


being 


sofa 


i= 


oe 


7x7 EL FH 


Compounds 8, 9, and 10 presumably have one o1 
the other of the above general structures, with the 
following substituents. 

8. 1,2:3,4-Di-O-isopropvlidene-L-arabinopyranose, 
m=; B=. 

9. 1,2:3,4-Di-O-isopropyvlidene-p-galactopyranose 
monoacetate, R=CH.OAc; R’=H. 

10. 2,3:4,5-Di-O-isopropylidene - p -fructopyranos: 
11 -C’-(hydroxymethyl) - 1,2:3,4-di-O-isopropyl- 
idene-p-arabinopyranose}], R=—H; R’=CH.OH; 
and the molecule is the mirror image of that 
depicted 

| ,2:4,5-Di-0-isopropvlidene-p-fructopyranose (com- 
pound 11) differs from the four other derivatives in 
this group in having a spirane structure at C2, instead 





er of a fused-ring structure. The boat conformation 
Pare « Reference groups * on the carbon atoms of the pyranoid rina of the exo-hoat confo-mation (IT) of some 1.2-(meth 
orthoacelates 
( ( ( ("4 ( 
(rrouy ] | , | G | 
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shown seems plausible; it somewhat resembles the 
ero form of the methyl orthoacetates (compounds | 
to 6), but differs in the position of the fused ring 





M1. 1.2:4,5-Di-O-isopropylidene-p-fructopyranose 


and in having the spirane structure. It may be 
noted that, whereas the pyranose compounds 1 to 
10 cannot exist in anomeric forms, a and 8 anomers 
of compound 11 should be capable of existence. 


3.3. O-Isopropylidene-Furanose and O-Isopropy]l- 
idene-y-lactone Derivatives 


Compounds 12 to 27 have either one or two dioxo- 
lane rings fused onto a furanose or a y-lactone strue- 
ture, or both Because these compounds have 
5-membered rings only, the molecules are more rigid 
than those having fused dioxolane-pyranose struc- 
tures. Each ring is substantially planar, and hence 
only one conformation seems probable. In com- 
pounds 12 through 21, carbon atoms 1 and 2 are 
shared by the rings, and in 22,5 23, and 24, earbon 
atoms 2 and 3 are shared. In each case, the fused 
rings are inclined at an angle of 120° to each other 

Compounds 12 to 19 have the following general 
structure 





12. |.2-0-lsopropylidene-p-xvlofuranose, 
R= HOCH R’ and R’’=H 


13. 1,2-0-lsopropylidene-p-glucofuranose, 


R—HOCH,-C -; R’ and R”’ =H. 
OH 

14. |.2-0-lsopropvlidene-L-idofuranose, 
OH 

R=HOCH,-C -; R’ and R’’ =H 
H 


15. Lithium 1,2-0-isopropyvlidene-p-glucuronate, 


R= LiO.C-C— : R’ and R’’=H. 
OH 
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16. Barium 1,2-0-isopropylidene-p-glucuronate, 


R= Ba,),0.C-C; R’ and R’’=H. 


OH 
17. Calcium 1,2-0-isopropyvlidene-L-iduronate, 
OH 
R=Ca,/,0,C-C—; R’ and R’’=H. 
H 


18. 1 ,2:5,5-Di-O-isopropylidene-p-xyv lofuranose, 


RR’’=H.CO——-C(CH,).; R’=H. 


19. 1,2:5.6-Di-O-isopropyvlidene-p-glucofuranose, 


R=H,C C R’ and R’—H 
() () 
C 
HC CHI 


The isopropylidene derivatives of b-glucuronolac- 
tone and L-iduronolactone (compounds 20 and 21 
have three fused, 5-membered rings; the furanose 
ring occupies a middle position, sharing Cl and C2 
with the dioxolane ring and C3 and C4 with the lac- 
tone ring. The dioxolane ring lies on one side of 
the plane of the middle ring at an angle of 120 
the lactone ring lies on the other side at an angle 
of 120°, but is not parallel to the dioxolane ring 
Compounds 20 and 21 have the following general 
structure 





20. 1,2-O0-lsopropylidene-p-glucofuranurono-6 3 -lac 
tone, R=H; R’=OH. 
21. 1,2-0-lsopropyvlidene-L-idofuranurono-6,3-laec 


tone, R=OH: R’—H. 


Compounds 22,° 23, and 24 have a dioxolane ring 
fused at C2 and C3 of the furanose ring and of the 
y-lactone ring, respectively. It may be noted that 
for compounds 22° and 23, @ and 8 anomers should 
be capable of existence. Compounds 22,° 25, and 24 
have the following general structure 





22. 2 .3-O-lsopropylidene-p-ly xofuranose,” 
R=HOCH,—; R’=H and R”’—OH 


or vice versa 








23. 2. 3:5, 6-Di-O-isopropvlidene-p-mannofuranose, 


R—H.C °° OCH, 
. 


rd 
H,( OCH; 


R’—H and R’’=QOH (or vice versa 


24. 2.3:5,.6-I i-O-isopropy lidene-p-mannono-1,4- 


lactone, R=H,C OCH,.; R’ 
Cc O-—O0. 
H.C OCH R’”’ 
The ketose derivatives 25, 26, and 27 have a 
dioxolane ring fused at C2 and C3: but if the 


compounds are regarded as substituted aldoses, these 
carbon atoms mav be classified as Cl and C2 of the 
furanose ring. Referring back to the skeletal 
structure V, these compounds have that structure; 
but, in each case, the molecule is the mirror image 
of that depicted. The following are the respective 
substituents 


25. 4:.¢ -O-lsopropy lidene -D ~thre O- pentulofuranose 


1-C -( hy droxy met hiv 1)-1 ,2-O0-isopropy lidene-p- 
thre o-tetrofuranose}, 


R and R’’ =H; R’=CH,OH. 


6 Di-O-isopropylidene - L-sorbofuranose 
L-¢ -(hvdroxymethy! -1,.2 3,5-di-O-isopropyli- 
dene-L-xvlofuranose}, 


26. 2,3 


RR’ R’ 


Potassium = 2,5: 4,6-di-O-isopropyvlidene-L-.rylo- 
hexulofuranonate, 


HoCOCNie,; CH.OH. 


2%. 


RR’ —H.COCMe,: R’— COOK 


3.4. Substances with a Fused Dioxolane Ring and a 
Fused m-Dioxane Ring 
Compounds 18, 26, and 27 have the following 


veneral structure 


In 26, R-CHLOH;: in COOK; and in 18, 
RH, and the molecule is the mirror image of that 
depicted. A m-dioxane ring (VI) is present in many 
sugar acetals; hence, determination of absorption 
which is characteristic of the structure is highh 
desirable | 


27, R 


ws 


3.5. Substances with ‘“Dangling’’ Dioxolane Rings 


The isopropylidene group attached to the oxygen 
atoms at C5 and C6 of compounds 19, 23, and 24 
forms a dioxolane ring that is free to swing about C4. 
Presumably, the position assumed by the group 
depends on the rest of the molecule, and the infrared 
spectra might be expected to provide information 
on this subject. Rings of this character, hereafter 
called “dangling,” should have properties quite 
different from those of the fused rings which are 
more rigid. 

| ,2:5,6-Di-O-isopropylidene-p-mannitol, compound 
28, has a dangling dioxolane ring attached by a 
carbon-carbon bond to each of the carbon atoms of 





28. 1,2:5,6-Di-O0-isopropyvlidene-p-mannitol 


ethylene glycol. Presumably, the spectrum of this 
substance is typical of spectra of compounds having 
dangling dioxolane rings, and it should be useful for 
assisting in the interpretation of the spectra of other 
compounds containing dangling dioxolane rings 
Although, for convenience, the ring can be con- 
sidered as capable of rotation (dangling), hydrogen 
bonding and steric factors may restrict the structure 
to a form in which each dioxolane ring lies in a 
definite sterie position. 


4. Discussion of the Spectra 


Comparison of the spectra of the compounds 
listed in table 1, on the of the structural 
features outlined in the preceding section, reveals 
that the absorption bands are not highly character- 
istic of the type of ring present, and that there are 
no readily distinguishable suitable for the 
assignment of ring structure or for the study of ring 
conformation of these compounds. On making a 
qualitative comparison of any one spectrum with 
another, it is seen, for example. that in a general 
way, the spectra of the mono-isopropylidene acetals 
resemble one another more closely than thes do those 
of the di-isopropylidene acetals, and vice versa, 
regardless of the particular sugar component or the 
size (or shape) of the sugar ring 

Although distinctive absorption bands are not 
apparent for pyranose or furanose, y-lactone, 1,3- 
dioxolane, and m-dioxane rings, surprisingly large 
differences are found in the spectra of substances that 
differ in what might seem to be minor details. These 
differences can be studied by direct comparison of 
the spectra for the closely related substances 


basis 


bands 


262 





4 





On account of the multiplicity of bands and the 
lack of highly distinctive ones, detection of absorp- 
tion characteristic for common structural features 
requires a statistical examination of the data. For 
this treatment, the compounds were classified accord- 
ing to their structural features, and the locations of 
the absorption bands were compiled. 

The method used for making comparisons is not 
conventional, but it seems suitable for the study of 
groups of related compounds. The distribution of 
the absorption bands for each group of related com- 
pounds was determined for the 2- to 15 u- (5,000 to 
667 cm range. First, the number of compounds 
showing one or more bands (or shoulders) in each 
wavelength increment of 0.1 « was counted, and the 
percentage of such compounds was calculated. This 
percentage was then plotted against the wavelength, 
to give the black regions of the plots in figure 1. 
Only one region (3.3 to 3.4 py) immediately revealed 
absorption by all of the compounds. Next, those 
groups of adjacent regions which seemed likely to 
encompass 100 percent of the compounds were noted 
and further examined; the results are given in the 
crosshatched regions of the plots in figure 1. 

Because frequency Is directly proportion: al to 
energy, use of wavenumbers in correlations of the 
kind sought is preferable to use of wavelength. On 
the other hand, the recordings (of infrared spectra) 
produced by the instruments we used are linear in 
wavelength, and readings from these charts must 
initially be made in terms of wavelength. Plots of 
the kind shown in figure | (which use 0.1-u inere- 
ments in wavelength) are compact, convenient, and 
reasonably satisfactory for use in analyzing a limited 
range. 

However, for such plots (linear in wavelength 
the bands appear crowded at low wavelengths (high 
frequencies) and become progressively less crowded 
as the wavelength increases (that is, as the frequeney 
The 0.1-u increments from 2.40 to 2.50 u, 
9.90 to 10.00 u, and 39.90 to 40.00 u correspond to 
decrements of 167, 10, and 0.6 em™', respectively. 


decreases 


From the practical standpoint, in making inter- 
comparisons, too much spreading out of bands is just 
as inconvenient as too much crowding of bands in 
such plots. Consequently, for the 10- to 40-u range, 
all wavelengths were converted to wavenumbers: 
and, in the analysis, decrements of 10 em=! in wave- 
number were employed. The resulting plots are 
not reproduced here, 

Ultimately, comparisons will be made by use of 
the classification system described previously |] 
using this statistical treatment. The present com- 
parisons are, however, based largely on the structural 
features characteristic of this relatively limited group 
of compounds. The compounds were studied in the 
groups outlined in table 3.) By making comparisons 
for groups of related compounds that differed in 

Various ways, the manner in which the absorption 
of the various compounds differed on the average 
was ascertained. Utilization of averages (instead 
of the individual absorptions) overcomes, at least in 
part, the “irregularities” found when the spectra of 


| 
| 
} 


| 
| 


individual compounds are compared. A few com 
parisons are given here; others will be deferred for 
inclusion in a study of a wide variety of carbohydrate 
derivatives. 


TABLE 3 tructural grou ps studied * 
(:roup Structural feature Compounds (numbers 
group 
l 0-Acety!]l (C=O 12.4to6,9 
2 Carboxyl, ionized (COO 15 to 17, 27 
3 1,3-Dioxolane ring : 1 to 28 
4 Hydrate 16, 17 
5 Hydroxyl group 3,7, 10 to 17, 19 to 23, 25, 2, 2s 
6 Lactone (C=O 2), 21, 24 
ry altro configuration 6 | 
‘S arabino configuration (ketopyra 1, 11 
nose 
4 arabino; galacto configuration (al 8: 6,7 
dopyranose 
10 Barium salt 16} 
11 Calcium salt 17 
12 m-Dioxane ring 18, 26, 2 
13 Furanose ring 12 to 23, 2 
14 threo; rylo; gluco configuration 25; 12, 18, 26, 27 13, 15, 1 
14, 20 
15 CGilycoside (8-pyranoside 5, tb 
Ith Hydroxyl group (glycosidi 22(?), 23 
17 Hydroxyl group (primary) 10, 12 to 14, ), 25, 26 
18 Hydroxyl group (secondary 3, 7, ll to 17, 19 te » 21, 22(? 
25, 28 
19 ido configuration 14, 17, 21 
yt) 1,2-O0-Isopropylidene acetal 7 to 21, 25 to 2s 
?] 2 .3-O-lsopropylidene wetal 22(?), 23, 24 
22 5,#-O-Lsopropy! idene acetal (“dan-| 19, 23, 24, 2s 
gle 
23 Lithium salt 5 
24 /yro; manno configuration 22; | to 4, 23, 24, 2S 
2 Methoxy! group or methy! ortho lto# 
wetate 
2h i}-C-Methyl (6-Deoxy) group 2.3 
27 Potassium salt 27 
2s Pyranose ring ltol 
2Y rue a-pyranose or 6-pyranose », 6, 11 
30 Pyranose ‘boat conformation? lto6, 11 
$1 Pyranose (chair conformation? 7 to 1 
32 Sugar ring l 23, 2 - 
$34 talo configuration j 
34 hree fused ring Sto 10, 18, 20, 21, 2 
* The bands found for the various compounds w he considered further in 


publications dealing with specific structural feature 


4.1. Absorption Bands Possibly Attributable to the 
Dioxolane Ring 


In 1949, Boekelheide and co-workers [7| reported 
infrared spectra for the 2-amyl-2-methyl and 2-hexyl 
derivatives of 1,3-dioxolane-4-methanol. In 1950, 
Kuhn [8] recorded the spectra (in the range 8 to 15 4 
of 1,2-0-isopropylidene-p-glucofuranose and its 5,6- 
0-isopropvlidene acetal (along with those of 77 other 
carbohvdrate derivatives), but made no attempt to 
analyze these particular spectra. In the following 
vear, Smith and Anderson [9] made a study of the 
infrared spectra of several 1,3-dioxolane derivatives 
110), and decided that it is “likely that the presence 
of a peak at 9.6 w indicates the presence of a dioxolane 
ring.”’ However, it remained for Bergmann and 
Pinchas [11] to conduct a detailed study of the 
infrared spectra of 18 compounds having the 1,3-ci- 
oxolane ring, and to make definite correlations; their 
data have also been interpreted by Barker, Bourne, 
and Whiffen [12]. Tsechamler and Leutner [13] later 
observed that compounds containing the 1,3-diox- 
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TABLE 4 thsorption bands 


Wavenumber 


Assigned by 


Smith and Anderson {1 
Bergmann & Pinchas tables |2] 
Bergmann & Pinchas text [2] 
Barker and co-workers [3] 
Pschamler and Leutner [4] 


formerly assiqned to the 


1190 to 1151 
1190 to 1LLSs 
1173 to 1151 


1 ,3-diorolane ring 


em of spectral range of absorption bands 


B Cc D> 


1042 
1053 to 103s 
1056 to 1038 


1105 to 1077 
1098 to 1068 


1151 to 1123 
1143 to 1124 


References 


1. L. I. Smith and R. Hl. Anderson, J. Org. Chem, 16, 972 (1951 

2 E. D. Bergmann and 8. Pinchas, Ree. trav. chim. 71, 161 (1952 

3. S. A. Barker, E. J. Bourne, and D. IH. Whiffen, Methods of biochemical analysis 
4. Ht. Tschamler and R. Leutner, Monatsh. 83, 1502 (1952 


olane ring generally show four recognizable absorp- 
tion bands. “Table 4 gives a résumé of the assign- 
ments of the prior workers. 

Our results are in fair agreement with the figures in 
the tables of Bergmann and Pinchas [11]: (a) 
Twelve of the compounds in the present study fail 
to show a band (B) in the LL51- to 1123-em~! 
(S.80+0.11 uw) region. However, Bergmann and 
Pinchas found that the exact location of this band 
depends more upon the alcoholic than upon the 
carbonylic residue of the acetal; in possible agree- 
ment with this conclusion, a band (B’) is observed at 
1126 to 1104 em™! (8.97 + 0.09 uw) for our compounds. 
(b) Two of our compounds show no band (C) in the 
1105- to 1LO077-cm~' (9.17 + 0.12 w) region. (ec) Eight 
of our compounds fail to show a band (D) in the 
1053- to 1038-em~! (9.57 + 0.07 uw) region. 

However, all of our compounds show one or more 


bands or shoulders in the (A) 1ISI- to 1153-em~'! 
(8.57+0.10 pu B’) 1126- to 1104-em™! (8.97 

0.09 w), and (€') 1093- to 1070-em~! (9.25+0.10 pw) 
regions. Except for the spectra of  1,2-0-isopro- 


pvlidene- and 1,2:3,5-di-O-isopropylidene-p-xvlofur- 
anose (compounds 12 and 18), all the spectra show at 
least one strong band (D’) at 1055 to 1029 em 
(9.60+0.12 uw for the two exceptions, there is a 
strong band at 1015 em™! (9.85 w). This character- 
istic of the p-xvlose derivatives is under further 
study; the spectra of the anomers and glycosides of 
p-xvlose are known to be exceptional [12]. 
four bands (A, B’, C’, and D’) seem to be the most 
important, but they are not particularly striking 
inasmuch as they appear in a rather crowded part 
of the spectrum. 

Each substance has an individualistic spectrum 
useful for identification purposes, but extensive gen- 
eralization requires great caution. Since many un- 
substituted and substituted sugars (and other organic 
compounds) also show absorption in 2, 3, or even all 
tof the above regions, allocation and Interpretation 
of these bands will be reserved for a more compre- 
hensive study (on a much larger group of sugar 
derivatives) to be conducted in the same manner. 

Many comparisons of limited scope can be made. 
For example, the spectrum of 1,2:4,5-di-O-isopropvl- 
idene-p-fructose (compound 11), a spiro—-dioxolane 
derivative, differs radically from that of its isomer 


These 


1151 to 11382 1105 to 1077 1053 to 103s 
doublet between doublet between 
1160 and 1120 1110 and LO) 
. vol. 3, p. 219 (Interscience Publishers, Ine., New York, N Y., 105 


(compound 10) and displavs 17 of the 19 bands 


recorded [14], in the range 2941 to 763 em for 
evclohexane-1 ,2-di(spiro-2’-1',3’-dioxolane). In ad- 


dition to these, it shows 26 other bands in the range 
3472 to 726 em". 


4.2. Other Absorption Bands 


The 28 compounds can be divided into three 
groups: the methyl orthoacetates, the mono-iso- 
propvlidene acetals, and the di-isopropvlidene ac- 
etals. Bands shown by all the members of each of 
these groups are given in the crosshatched regions 
of figure 1 (B, C, and D) and in table 5. It is noted 
that, in certain regions, more than one interpretation 
of the data is possible. 

Bands shown by all of the compounds in table | 
are shown in the crosshatched regions of figure 1 (A 
and in table 5. Selection of one of the alternative 
ranges exhibited by a subgroup could often be made 
by intercomparing the ranges for all three subgroups. 
(Even so, more than one interpretation of ranges for 
band locations was found possible in certain regions 
of the spectrum.) 

In addition to those mentioned in the preceding 
the following are of the bands 
shown by a// of the compounds, with possible assign- 
ments given in parentheses: 3012 to 2967 em~!, Le 
3.32 to 3.37 uw (CH stretching); 1468 to 1439 em 
ie., OSL to 6.95 uw (or 1466 to 1437 em 6.82 


and 1333 to 1312 em~'. Le... 


section, some 


to 6.96 w), 7.90 to 7.62 wu 


(C—-H_ bending): 1391 to 1379 em~-'. ie... 7.19 to 
7.25 pw (or 1385 to 1372 em, Le., 7.22 to 7.29 yp), 
often as a doublet (possibly — CH, deformation, but 
these bands are shown by many unsubstituted 
sugars); 1272 to 1236 em™', Le. 7.86 to 8.09 pg 


(C—O), and 1244 to 1212 em™', i.e., 8.04 to 8.25 uw. or 
1232 to 1205 em™!, Le., 8.12 to 8.30 uw (C—O stretch 
ing); and 889 to 858 em~', Le., 11.25 to 11.66 uw, on 
877 to 850 em~', 1e., 11.40 to 11.76 w (1.3-dioxolane 
ring?). Other bands, which may possibly be corre 
latable with one or other of these features, were 
noted at 981.4 to 948.8 em~', 824.4 to 793.0 em 

and 540.5 to 507.6 em Some other possibilities 
noted encompass broader ranges of wavenumber and 
show overlapping. (They are mentioned in 
some of them should prove be significant, but 
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PABLE 5 
W ivelengt! ul 
Phe six methyl rhe eleven mono 
orthoacetates opropylidene acetals 
me, Te 2 ‘ 
§. 33 to 3.37 §. 32 to 3. 37 
or 3.34 to 3. 3Y or 3.34 to 40 
6.75 to 6.84 o1 
6.81 to 6.95 6.82 to 6. ¥2 
or 6.85 to 6. 46 6.82 to 6. WZ 
6.95 to 6.47 
1.19 to 7.23 20 te 2 
7.26 to 7.20 or 7. 22 to 7.27 
7. 52 to 7. 5S 7. 0 to 7. 62 
64 to 7. 88 
2 te W7 ol 
7.43 to 8.09 7.86 tos. 04 
Wi tos. OY or 7.41 to 8. 12 
or 8.04 to 8.15 rs.04 to 8, 22 
& lO to S. 26 
ors. 12 to 8. 30 or 8. 12 to 8, 30 
& S44 to 8.64 x 4 mf 
8.72 to 8.79 
s Ys tov ys ( i 
4W 15 to 9. 30 UY 2WtoY 
4 34t09.4 G2 1390 
4 35 to 4. 4s 
4. 62 to 9.79 +61 to OS 
uH2to &7Y or 4.72 to ¥St 
10. 10 to 10, 21 or @ &5 to » 
rlOWto WS we. 321 lf { 
10.74 to O80 
10.93 to LO 
11.08 to 11, 1S lll to bb 3e 
11.40 to 11.61 or 11.26 toll 
11. 40 to 11.61 11. S2to lit 
12.25 to 12.34 
urther work will be needed in order to test this 
pot These ranges are: 664.5 to 612.4. em 
580.4 to 524.9 em 555.6 to 518.1 em; 532.8 to 
193.8 com 530.8 to 487.1 em 6.7 to 38767 
em 128.5 to 365.6 em 118.4 to 363.0 em 


$110.2 to 349.4 em and 370.1 to 319.0 em 
The following absorption bands were shown only 
having the structural 


features mentioned (groups 1, 2, and 4 to 6 of table 
the possible 


3): the and 

assignments are given in parentheses. The acetates 
1748 to 1739 em™ (5.72 to 5.75 wp: 
and at 607.9 to 598.1 em™! 
Those comM- 


(r 


by compounds (in this study) 


wavelengths in microns 
showed a band at 
acetate C—O stretching 
(16.45 to 16.72 uw: acetate deformation 
pounds having a hydroxyl group and an ionized 
carboxy! vroup showed bands at 3145 to 3125 em 

OH modified by C—O), at 2786 to 
OH modified by C—O 


(3.18 to 3.20 pw: 


2752 ¢m 3.99 to 35.606 uw: 


and at 1637 to 1600 em (6.11 to 6.25 uw: COO 
The hydrates showed a strong shoulder at 1645 em 


6.08 w; water of ervstallization) and a band at 706.2 
to 703.7 em All of the compounds having one or 
more free hydroxyl groups showed a band at 3472 
to 3279 em (2. SS to 3.05 UL: associated alchoholic 

QO —-H, stretching frequenes Compounds 10 and 
28 showed a band at 3289 to 3279 em™! (3.04 to3.05 u: 
H-bonding The lactones showed a band at 1789 
to 1764 em 9.909 © stretching 
queney 

All of the spectra were also studied in regard to 
the 28 other features listed alphabetically as groups 
7 to 34 of table 3 Those compounds having a 
m-dioxane rine showed a band at 310.3 to 307.7 


to 5.67 uw: C Ire- 


Bands shown by all of the 1.38-dioxrolane derivatives and by all the members of each of the subg OW ps 


pectral range of absorption bands 


rhe eleven di-isopro All of the 28 acetals 
pylidene acetals 
82 to 3. ob 3. 32 to 3. 37 
6.83 to 6. 87 6.81 to 6. 05 
6. 83 to 6.87 or 6. 82 to 6. Ft 
7. 20 to 7. 25 7.19 to 7. 25 
or 7. 22 to 7. 28 or 7. 22 to 7. 29 


7. 30 to 7. 59 7. WH to 7. 62 


7. 6S to 7.83 


7.88 to 8.02 7.86 to 8.04 
& 15 to 8. QA or 8.04 to 8. 25 
8.15 to 8. 25 - 

S.15 to 8. 25 or 8. 12 to 8. 30 
or 8. 20 to &. 31 

SS. 56 to 8. 67 s. 47 to 8.67 


s Sp TOS. US OT 


SSS TOU 04 SSS TO 4. Ub 
9, 22 to 9. 35 Y 15to¥ 
9, 23 to 9. 37 or 9.23 to &. 45 


4 SO to & S85 
or 4. Ad » USS ora 


4. AO to SS 
55 to 4 SS 


or ¥ S82 To 4 US 
4 SS to 10.06 or 4. S5-to 10. 22 


10. 283 to 10. 43 or 10.19 to 10.44 
10. 7S to 11.04 


lL. 


Ll. vt 


5, to 11. 66 11. 25 to 
r 11.51 to 11.76 or 11, 40 to 


cm The 3-glycopyranosides showed bands at 
967.1 to 966.2 em~', 811.0 to 810.4 em 646.0 to 
644.8 ¢m and 475.1 to 470.6 em Those com- 


pounds having a free hydroxyl group at carbon 
atom | showed a band at 893.7 to 892.1 em”! and 
at 683.1 to 680.7 em The isopropylidene acetals 
showed a band at 859.1 to 833.3 em~'. Those com- 
pounds having a sugar ring (including compound 24 


showed one or more bands at 914.1 to S8S4.2 em 


(Other possible interpretations are: 900.9 to S7S8.7 
ecm” !, or 900.1 to 877.2 em‘, or 896.9 to 867.3 em 
or S92.1 to S6L.5 em They also showed bands 
at 653.6 to 611.6 em and at 591.7 to 544 em 


Other possible ranges (which, because of the broadet 
range or the overlapping, may not be significant 
are: 639.0 to 578.0 em 584.8 to 540.5 em and 
580.4 to 533.3 em Further work may throw light 
on the significance of these bands 

From a study of the infrared absorption spectra of 
sugars and their derivatives capable of existence in 
anomeric forms (namely, of compounds in which the 
hvdroxvi group of the hemiacetal carbon atom is, 
for example, free, acetylated, or involved in gh 
cosidic union), Barker and co-workers |15, 16} were 
able to recognize certain absorption bands which 


o}) : 


were considered to be characteristic of the @ and 2 
respectively, of the pvyranoid forms. In 
pvranoid compounds 
nnomers are 


nnomers, 
the present study, the only 
capable of existence as a end B 
and 11 (which are, respectively, a 
b-p-galactoside, and a p-fructo- 
, and they show bands fitting the 


com- 
pounds 5, 6, 
B-p-glucoside, i 
pyranose derivative 
assignments called types 1, 2, and 3 by Barker and 
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co-workers. The other pyranoid derivatives, which 
cannot exist in anomeric forms, show many of the 
bands in question; but, as might be expected, these 
do not fit the structural correlations. 

Similarly, of the furanoid derivatives, only com- 
pound 23 (and possibly * 22) is capable of existence 
in anomeric forms, and it shows bands fitting the 
tentative assignments (for furanoid sugars) called 
types A, C, and D by Barker and Stephens [17]. 
The other compounds show many of these bands, 
but they do not fit the structural correlations. 


5. Experimental Procedures 


5.1. Preparation and Purification of the Compounds 


The compounds listed in table 1 were prepared 
by the methods given in the references cited. For 
compound 4, the few milligrams of the material 
remaining from the original preparation of Pigman 
and Isbell [3| had a melting point somewhat low, 
but it was apparently in good condition; because of 
the lack of material for recrystallization, its infrared 
absorption spectrum (2 to 15 yu) was determined 
directly, and this showed the presence of a small 
proportion of free hydroxyl group (band at 3484 
em'; 2.87 uw). For compounds 21 and 25, the 
melting points were satisfactory, but insufficient 
material remained for checking their optical rota- 
tions. With these exceptions, each substance was 
recrystallized or otherwise purified until further 
treatment caused no change in its melting point or 
optical rotation. In most cases, these values then 
agreed with those found in the literature but where 
a higher melting point is now given, it is believed 
that the new values reported in table 1 are correct. 

As regards the preparation of compound 12, by 
the acid hvdroly sis of compound IS, a survey of the 
literature showed that use of 0.16 percent hydro- 
chloric aeid {18, 19], 0.2 percent hydrochloric acid 
20, 21), 0.04 N sulfuric acid [22], and 0.2 N sulfurie 
acid [23| had been recommended; the last concentra- 
tion cited may be a misprint. Accordingly, a solu- 
tion of compound 18 (10 g) in 250 ml of 0.05 N 
sulfuric acid was kept at room temperature for 2 


hours. The acid was then neutralized by adding 
barium carbonate (5 ¢) and shaking until neutral 


to Congo Red. The suspension was filtered through 
a laver of Celite plus barium carbonate, and the 
clear, colorless filtrate was added to barium carbonate 
(1 ¢) and evaporated to dryness under diminished 
pressure at 30° (bath temp. The mixture of sirup 
plus solid was extracted with acetone and filtered, 
and the filtrate was dried with anhydrous sodium 
sulfate, filtered, and evaporated to dryness under 
diminished pressure. The resulting sirup was dis- 
tilled (bp 112° to 114° © at 0.06 mm), affording an 
almost quantitative vield of a colorless, viscous 
sirup) which completely ervstallized spontaneously. 
Although previous workers [18, 23] have stated that 
the compound cannot be recrystallized from organic 
solvents, it may readily be recrystallized from 
ether (2 vols.) plus heptane (0.2 vol.) on standing 
in the freezing compartment of the refrigerator. 


5.2. Preparation of the Pellets 


Samples for spectrophotometric study were pre- 


pared in the solid phase, as pe ‘lets of the ¢ rystalline 


compound suspended in an alkali-metal halide. 
The pellets were 9.5 mm in diameter and weighed 
approximately 100 mg; the weight of sample per 100 
mg of potassium chloride (for the 2- to 15- wu region) 
or per 100 mg of potassium iodide (for the 15- to 
40-u region) is shown in table 6. The accuracy of 
we sighing was not high (40.1 mg per mg of sample ), 
but this was woe at 5 unimportant in view of 
other factors to be dise cussed. No spec ial pre- 
cautions, beyond working under an infrared lamp and 
storage in a vacuum ‘desice: ator, were taken for 
drying the powder, which was hand-ground in a 
Mullite mortar. The pellets were pressed in a simple, 
nonevacuable, three-piece die, consisting of a eylin- 
drical body and two plungers (9.5 mm in diameter) 

similar to one described by Anderson and Woodall 


[24]. A machine-shop vise was used for applying 
pressure (estimated to be 20 tons/in2 for pellets of 
this diameter). Satisfactory pellets were obtained 


without evacuation. It should be noted that there 
may have been a slight variation in pellet thickness 
from pellet to pellet. 


TABLE 6 Conditions for recording the spectra 


Com Molecular Concentration for 2-to Concentration for 15 
pound weight l5-u range to 40-4 range 
mg 100 m@ mole/ka of mq 100 m@q mole ka of 
AC! mirture Al nirture 
l s62. 35 o.4 O01; yd 0.05, 
2 404, 30 } Ol 2 a 
3 220. 23 4 Ols 
4 362. 35 4 Ol 
5 O50. OO 4 OO) 2 osu 
i O50. 00 $ O06 ” O30 
7 220, 23 5 02 2 l 
~ 230, 26 } ol » Os 
u 302. 33 1 Ol r i) 
10 20, 24 } vl 2 “7 
11 20). 2 ; 01 e i 
12 140. 20 ) OS i 
220. 23 } Ols 2 Us 
22). 23 4 Ols 2 OS 
Fa) 240,14 4 01 2 Os 
It) * S10, 8U 4 Ol ; (Wy 
7 271. 26 i * Ol - UF 
Is 230), Jt l 4 pi Os 
lu 2), 24 u.4 ol 2 i, 
yal 216. 20 } Ol« » au 
yA | 216, 20 4 Ole 2 ay 
22 190), 20) 3 17s } yal) 
23 20, 20 “4 ol J OF 
24 JAS, 2 l OSs« 2 “ 
25 140. 20 0.8 4 j A“) 
2h 20). 24 i Ol e i 
27 312. 37 l 03 oe) 
28 22. 31 } Ol r 
Equivalent w h 


5.3. Measurement of Infrared Absorption 


The spectrograms are shown in figures 2 and 3 
Those in the 2- to 15-4 range for compounds 7, 12, 
IS, 22, 24, 25, and 27 were recorded with a Beckman 
Model IR4 (double-beam) spectrophotometer equip- 
ped with prisms of sodium chloride; they have delib- 
erately been reproduced somewhat differently from 
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the rest of the spectra All the other spectrograms 
were recorded with a Perkin-Elmer Model 21 (double- 
beam spectrophotometet equipped with a prism of 
sodium chloride (for the 2- to 15-u range) and of 
cesium bromide (for the 15- to 40-u range The 
use of these two prisms was found to be adequate for 
obtaining survey or identification spectra encompass- 


ing the entire 2- to fO0-u range, because we have 
found that the resolution obtamable with a cesium 


bromide prism is not much inferior to that given by 
a potassium bromide prism in the 15- to 20-~ range 
and is about its equal in the 20- to 25-y range. By 
replacing the usual, stray-radiation filter-turret with 
containing retractable filters (Reststrahlen 
mirrors) of lithium fluoride (for the 15- to 30-u region 
and of calcium fluoride (for the 30- to 1O-u region 

the Stray radiation (encountered in the 15- to 40-u 
region When a prism of cesium bromide is used) was 
virtually eliminated for the 15- to 34-u range, and 
reduced to one percent at 


one 


$6.54 andto two percent 
at 38 uw. No pellet blank was inserted and, indeed, 
this is not at present practical for the far infrared, 


because insufficient energy is available for that 
revion 

In the presence of strav radiation, it is not in general 
possibl oO auchtie complete compensation for substances 
such as water vapor ommon to the reference and sampl 


ve iis OF the spectrophot 


The 


reference 


recorade d 
re fle ection 
ibsorption Dv traces of moisture 


hiealice Con- 


oOmetel! 
} 


compensator nm the 


spectra were 
without beam for 
rseatter bv the pellet. or for 


sugar derivative or in the alkali-metal 
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observed il 2 


in the 2- to l5-u range 
attributable tO water 
Attention is drawn to these regions by me 

Similarly, 


wate! 


sequently ibsorp- 
and 6 | MM. 
tns of dashed lines 
UNCOML pe sate d absorption by atmosplhe ric 
was evident in the far-infrared curves, and the 
spectral region bevond about 344 was doubtful in the 
that weak bands might not be 
imposed 


tion Wiis 
SOTTLE 
Vapol 
Se LSE 
observed beenus ol thre 
bands due lence, 
spectrograms dushed 
shown in fig. 4 


sUuper- 
uncompensated to wate! these 

v drawn on the 
inserted after 
obtained for 


Nujol. These 


should not be 


regions are with lines, 
consideration of the 
potassium iodide, moist potassium iodide, and 
dashed lines are merely pre 


interpreted quantitatively 


CUrVes 


mtionary they 


In the 2- to 15-u range, the accuracy of the wave 
length readings was better than 0.02 uw: in the 15- 
to 40 w range, the accuracy was within =—0.05 u 

Because of the possibility of interaction of the 
Various compounds with the pelleting halide under 
high pressure, the spectra of typical compounds were 
also recorded in a Nujol mull (requiring no pressure) 
for comparison. | 
visual (development of color) and speetrographic 
evidence of interaction with potassium iodide and, 
to a smaller extent, with potassium bromide, but not 
with potassium chloride; however, this behavior was 
not observed on repetition with fresh samples of the 
halides 

The spectra (2 to 15 w) given for compound 18 
illustrate how exhibited Iyy a compound 
suspended 11) potassium chloride can be obscured (by 
the bands of Nujol) when a Nujol mull is employed. 
Another reason for preferring the alkali-halide pellets 
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is that some of the compounds studied are not suffi- 
ciently soluble in such solvents as dioxane, carbon 
disulfide, chloroform, or tetrachloride to 
afford solutions concentrated enough for recording 
ol satisfactory infrared absorption spectra 


carbon 


In table 6 are shown the calculated concentrations, 
in mole per kilogram of alkali-metal halide, of the 
These varied rather widely because 
of (a) the range in molecular weight of the com- 
pounds and (b) the different concentrations (in 
milligram per milligram of halide) which were found 
for adequate recording of the spectra 


Various pellets 


hecessary 
\s a consequence, COMparisons of intensity of ab- 
sorption, from one compound to another, can only 
be true and quantitative where the concentration 
(in mole per kilogram) Intercompari- 
sons made for the other compounds are relative and 
qualitative as regards the intensity of the 
although the positions of the bands (in microns) 
Trhany properly be compared. 


is the same. 


bands, 


Insufficient material was available for recording 
the spectra of compounds 3 and 4 in the 15- to 40-y 


revion 
5.4. Spectra Measured Under Different Conditions 


The infrared ervstalline 


materials show more bands than the spectra of the 


absorpt On spect ra of 


Sine compounds It) solution This Is because, for 


the solid, the following factors are operative l \ 
removal of degeneracies by perturbations from lattice 
forees, (2) an intensification of overtones and com- 


bination Vibrations, and (3 occurrence of combina- 
tion vibrations between molecular and lattice modes 
in the solid | 

The the 2- lo-u 
pounds land 2in carbon tetrachloride and in dioxane 


had been published previously [25]; those given in 


spectra (in to region) of com- 


pellets Of potassium chloride for the same region of 
wavelength are shown in figure 2 \ comparison of 
the relative complexity of the spectra for the solu- 
tions and pellets of compounds | and 2 (as well as 
of many other sugar derivatives vet to be published 


shows the above-mentioned effect 


The authors express their gratitude to Harriet L 
Frush, J. D. Moyer, and R. Schaffer for providing 
samples of certain of the compounds used in this 
study; to Alex Cohen for technical assistance; and 
to J.J. Comeford and F. P. Czech for recording some 
of the infrared absorption spectra. 
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Earth Currents Near a Top-Loaded Monopole Antenna With 
Special Regard to Electrically Small L- and T-Antennas 


H. Lottrup Knudsen 


An investigation has been made of the ground Currents near a top-loaded monopole with 
nonazimuthal symmetry Formulas have been developed for the surface current density 
produced by an inclined, straight wire over a horizontal ground plane for an arbitrary current 
distribution on the antenna Working formulas have been developed and numerical caleula- 
tions of the surface current density on the ground plane have been carried out for the case of 
a small antenna with a linear current distribution These results have been used for the 

ileulation of the contribution to the surface current density due to the top loading in the 
case of an L-antenna and in the case of a T-antenna both the absolute value of 
the surface current density arising from the top loading and the relative value of its ¢-com 


In ( ach Cust 


ponent have been plotted, us it 
lnportant 


Thay 
cumstances mav be 
radial ground wires 


1. Introduction 


Low-frequency and very-low-frequency antennes 
often consist of a vertical wire with a top loading 
of inclined, or in a special cese, horizontal, wires 
radiating from the top of the vertical member as 
shown in figure la 
the ground around the antenna, caused by the cur- 
rents induced in the earth, a system of ground wires 
is generally used. As the ground currents have 
essentially a radial direction, it is a common practice 
to let the ground wires point radially outwards from 
the base of the antenna. 

In the far zone field of the antenna the main con- 
tribution to the field comes from the vertical wire 
If the top loading is horizontal, the ratio of the con- 


In order to reduce the losses in 
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oaded antenna by issu 


FIGURE |] 1, Top 


hution on this antenna 


in caleulating 


be expected that this component under certain cir 


case of a svstem of 


he gro ind losses in the 


tribution to the field coming from the horizontal 
wires to the contribution coming from the vertical 
wire tends toward zero, as the distance from the 
antenna goes toward infinitv. Also in calculating the 
ground currents and the ground losses, i.e., in ealeu- 
lating the field at the ground near the antenna, it has 
generally been assumed without proof, that the con- 
tribution from the top loading is insignificant Re- 
cently Wait! has investigated this question. In his 
investigation he assumes that the top loading may be 
in the special case of a hort- 
zontal top loading, by a disk; i.e., he assumes perfect 
rotational symmetry In his Wait 
arrives at the result that under usually existing cir- 
cumstances the top loading contributes only to a 
small extent to the currents-and losses in the ground 

However, as has been pointed oul by Wait, when 
the number of wires in the top loading is very small, 
sav one, as in the case of the L-antenna, or two, as 
in the case of the T-antenna, the top loading may 
furnish a circumferential component of the ground 
current, which, although small, may contribute con- 
siderably to the ground losses, as it is perpendiculas 
to the radial ground wires. It therefore, 
justified to carry out an investigation of the ground 
current excited by a top loading consisting of only a 
few wires and to compare it with the contribution 
from the vertical wire This is the object of the 
present. paper. Whereas formulas for the ground 
current have been worked out for the general case, 
working formulas and numerical calculations are 
given only for the case where the antenna is small 


replaced hy a cone, OF, 


investigation, 


Secs, 


compared to the wavelength. 

The theory has been developed under the dssulnp- 
tion that the current distribution on the antenna 
including that on top loading is known. For an- 
tennas that are short compared to the wavelength, 
we may make the following statements: If the num- 
ber of wires in the top loading is large, the current 
distribution will probably be very near quadratic 
Wait in his investigation of a 
the other hand, if the 


as was assumed by 
disk-loaded antenna. On 
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number of top-loading wires is small, say one or two, 


as we shall assume here, the current distribution 
will probably be linear, to a good approximation 
This is shown tn figure Ib, where A denotes the 


length of the vertical wire, a, the length of the radial 
wires, (these being assumed to be of equal length 

and J, the antenna current. A part of the following 
theorv has been derived for an arbitrary current dis- 
tribution on the antenna. However, in those parts 
of the investigation w here a specific choice of the ecur- 
vo. in the nu- 
distribution 


rent distribution has been made, as e 
merical calculations, the linear current 


| 
has Hpeenh ¢ nose 


List of Principal Symbols 


a—length of radial wires in top loading, 


A-vector potential of current on inclined 
Wire, 
| coefherent in power 


current distribution funetion ¢ 


SCTICS 


CNPansion lor 


/; height over ground plane of one end of 
inclined wire, 


0 h 

IB sinh a Cos (od re} COS a, 
( c 

( leneth of inclined wire, 
D h 

C=(-)+(-): 


( c 
arbitrary 
function for 


7 reference distance 


distribution inclined 


f(s current 
wire in terms of s 
do current distribution function for imelined 
wire in terms of o, 
/ length of vertical member of antenna 
/ current on inclined wire, 
reference current (arbitrary 
/’ current in one top-loading wire at 
nection point with vertical wire, 


/ current in vertical member at ground, 


Cort 


/ eurrent in vertical member nut connector 
: point with top loading, 
normalized surface current density on 
eround plane ] A 
normalized surface current density On 


eround plane excited I one horizontal 

Wire, 
normalized 

eround 


surface current density on 


plane excited by horizontal 
] -antenha, 


current 


members of 
normalized surface 
ground plane excited — by 


density on 
horizontal 


members of antenna with four top- 
n loading wires, 
normalized surface current density on 


eround plane excited by vertical mem- 
ber of antenna with top loading, 
fine 


integral ds. 
s / 


oJ = the 


Av surface current density on ground plane, 
nv —number of wires in top loading, 





N=( I, 


p—moment of Hertz dipole, 


p the integral ( ( z ) — ike ( - ) | de. 
P the integral ( ( 3) ike( 5 ) | a)e*"*de, 


r—distance from point of antenna (£, n, ¢) to 


field pont (p, 9, , 
P— distance from point of antenna &, nn, ¢) to 
point of ground plane (p, ¢, 0 
N coordinate along inclined wire, 
S unit vector pointing In positive direction 
of inclined wire, 
unit vectors pointing in the two directions 
T-antenna, 
inclined wire 


of horizontal members of 
S—pomt of penetration of 
through ground plane, 


ret unit vector in z-direction ete,, 
\ R? 
eat 
\ ( oF 1.7 


a angle between downward vertical direction 
negative c-axis) and positive direction 
S$) of inclined wire, 
3—azimuth of inclined wire, 
anv one of the coordinates p, O Or ©, 
evlindrical coordinates of field point, 
0 normalized coordinate s ¢ along iIncelined 
wire, 
En, € rectangular coordinates of pont of inclined 


wire 


2. Earth Currents Excited by Inclined Wire 


In this section the earth current excited by the 
current in an inclined straight wire will be calculated, 
with specail reference to the case where the leneth 
of the wire, as well as its height over the ground 
plane, is small compared to the wavelength 

The vrounad plane is taken as the (er, y plane in 
A wire of the length 


an (wr, ¥, coordinate svstem 

e with one of its ends situated on the z-axis at the 
herht h over the cr plane forms the angle a 
with this axis as shown in figure 2 A coordinate 


s along the wire is introduced, s=0 corresponding to 
the point of the z-axis. The wire is assumed. to 
carry a current J(s)—J,f(s) where J, is an arbi 


trarily chosen reference current. In this section no 
particular choice of J, is made 

The normalized coordinate sc and the function 
W\o (sie f(s) are introduced The antenna ay 
be deseribed in parametric form in the rectangular 


coordinate svstem (7, y, 2) as follows 


< S Sill @ COs OO, 
n S SIll @ Sill oO, (‘)}=~s ( 
C—b—s cos a, | 
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FIGURE 2 / ned e above ground plane 


With the time factor the reetangular coordinates 


of the vectol pote ntial \ are expressed by 


ul 
| sina cos B J, 


| sinasin B J, 
lr 
| COS a a 2 
where 
/ \ ‘ 
e ds , 
: j 
with 
si! ¥r COS O pop COS 0 SSilnasin ob 
p sll @ - h S COS @ : } 
where p. @ denote the evlindrical coordinates 


of the field pot The coordinates of the 


potential Ain the evlindrical coordinate system are 


vector 


ul 


A | COS @O el, sin 78) Sill a Cos (od Db m i 
ul 
A 1. sin @-.1, cos @ 4 sin a sin (@—8) J. 
ul 
A cos a J 5 


The Inavgnetie field intensity I/ may now be obtained 
from 


“1 


285 


Until now, the effect of the ground plane has not 
been considered. Assuming this to have infinite 
conductivity, this may be done by adding the field 
of the image. Since we are interested only in the 
tangential magnetic field at the ground plane, this is 
twice the tangential component of the magnetic 
field calculated above 
For the current density A in the ground plane, 
then, 

A==Z /] 


(0) 


[It is convenient to define a dimensionless 


, hormalized 
surface current density 


d— a 
yh 


where dis an arbitrarily chosen reference distance. 


With the above definition of 7 it is found that 


Defining 


SSI a@ cos 3 pcos © ‘sll asinp 


YS COS a by 


= P 
c ( o SIlLa@ COs bm COS ¢) +( o sin @ sin re 
e 
0 h 
sin @ ) ; ( G6 COS an ) ’ 
‘ \ 


we finally obtain after some manipulation 


p sill © 


| d h ) 
» | - Sil @ cos (@— PB) 4 cos a |P, Lla 
gu Cc Cc Ya 

| f h ° 
= SIN @ SIN(® 8) P. llb 


The first term in the square bracket in the expression 
for P corresponds to the induction field, whereas the 
last term corresponds to the radiation field 

It is easily verified, from the above expressions for 
the components of jas well as by physical considera- 
tions, that jp 8 an even function of @ 
is an odd function of ¢—<, 

As a partial check of the formulas developed here, 
these formulas a be applied to a short vertical 
wire with a constant current / 


3, whereas }, 


as considered by 








Wait see Tootnote | This Tian \ be done, for example, 
by setting 
h () 
q o 
‘ f ow 4 12 
Then 
5 / / p Le p 4 | ik: 
/ COS a, a : ds 
: d Dr ¢* lr ‘ E P? 
13 
WY ()} 
where 
R lp- y¢ 


Phis expression checks with formula (7) in the above 
mentioned paper by Wait (except for the sign; a dif- 
ferent sign convention seems to have been used by 
Wait 
Whereas it 


It) order to 


to evaluate the integral 7? 
normalized current 


Is hecessary 
the 

density 7, the direction of 7 may be obtained without 
earrving out anv integration. Thi this is 
that all the current elements of the wire give con- 
tributions to 7 pointing in the same direction. The 
direction of the current is expressed by the angle 


obtain surface 


reason for 


Which it forms with the p-axis, as shown in figure 3 


It is found that 


h Slli@ sill (@ Oo 
tain ; 
YO Sill a@ COS (OD Db 


lt \W ill be shown threat this eXpression exactly expresses 
that the surface current density y,yat every pont of 
the ground radially outward from 


lane OES 
| | 


the pot S in which the inclined wire 
or its prolongation penetrates the ground plane 


or toward) 


Referring again to figure 3, it is found that the 
distance from the origin O to the point 5S is 6 tan a 
Introducing the perpendicular SM from S to the line 
from the origin O to the field point N, if the current 
at N points radially outward from 5, it follows that 
S\ b tana Sih (@ re] = 

tan . Lo : 
NM p—b tan a cos (¢—8 


u tan 





As this expression is equal to the expression for 1 
derived above, it is concluded that the surface 
current density points radially outward from = (or 
toward) the point where the linear antenna or its ‘ 
prolongation penetrates the ground plane. This | 
result, which was arrived at as a detail in a more | 
comprehensive investigation, may of course be de- | 
rived in a more straightforward manner. In figure 
fa is shown an antenna that penetrates the ground 
plane at S, and lines that are everywhere parallel 
to the surface current densitv. Thev are not field 
lines in a strict sense, as their mutual distance is 
not proportional to the magnitude of the surface 
current densitv. In the case where the antenna is 
parallel to the ground plane, i.e., Where the pont ol 
penetration is at infinity, the surface current density 


will be evervwhere parallel to the wire. This is 
illustrated in figure 4b 
y y 
§ ‘ 
* 
lane aT ¢ 
* ; FA 
eA + ~ -~ x — st -_ x 
Q b : 
FiGuRi j Line how lirectio t } “a 
density on the ground plane due toa near antenna 


; 


Returning to the ealeulation of the magnitude — of 
the surface current density, the expression 7? must 
be evaluated In the veneral case the Integral by 7 
which 7? is defined can apparently not be expressed 
by known funetions In the general case. therefor 
numerical or graphical methods must be applied to 
evaluate ?. This can be done in any specified case 
without a prohibitive amount of labor 

A case which has a great practical importance and 
in Which the analysis is greatly simplified is the on 
where both the height 6 and the length ¢ of the an- 
tenna are small compared to the wavelength or more 
exactly where 


hib+e | 
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A similar approximation was used by Wait (see 
footnote 1) in the case of a vertical antenna. For 
an inclined wire, as is considered here, this approxi- 
mation will, for the directions perpendicular to the 
wire, be just as good as for the vertical wire. How- 
ever, in the directions of the projection of the wire 
on the ground plane, the approximation will be worse 
than in the above mentioned case. Or otherwise ex- 
pressed, for the directions not perpendicular to the 
wire, ke must be smaller than for a vertical wire in 
order that we shall obtain the same approximation. 

The rest of this report will be concerned exclusively 
with this case. Under the above assumption, any 
case will be included in at least one of the following 
four situations 


l. hp 1, 

2. Neither kp l nor be p, 
s be » but not | hp, 

1. | hp 


When Aib+e) is not much smaller than 1, say 


k( be) 0.2, situation 2 may apply for some values 
of p. On the other hand, if #(b+e 1, situa- 
tions | and 3 may apply simultaneously for some 
values of p The above discussion is illustrated in 


figures Sa and Sb applying respectively to the situa- 
tion where k(b—e) is relatively large (although con- 
nequality Lib+e 1) and to the 
is extremely small, & h c 


forming to the 1 
situation where /&(b— ¢ 














k (bec)< 
Jeither k +c<«< 
p<< nor bec “op yf Ke Ss“ Kg 
< o-~ o ~ 
_-_-—- Ko 
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+ 
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o ~ i Kp KP 
a ole 
_-—-— Kp 
U 
k (b+ 
bh 
— 
RIGURE 5 Relat position of the ariou eqions of paran 
/ ; ; , 
; } } approcimations to he nieara ] 
ipp 


kp | 


pr | 


1] ih I ihkh T = “|s 0 ido 


~[.Ge) Lea" arden | (ip) store 


terms containing (ka)*, (kh)?, (kp)?, ete. (16 


In using the last approximation the radiation field 
term, as Well as the retardation, is neglected. 
Neither kp | nor be p 

It turns out that the last one of the approximate 
formulas worked out for the case of kp lis a 
special case of the approximate formula in the case of 
Probably, the last formula may there- 
applied with a fair approximation in the 


present cnuse, too 


h r¢ p. 
fore be 


O--Cc p but not 1 hp 
> ‘y S - ’ 4 ) S b 
I p l é SI @ COS (dO re +( ) é COS a 
p p Pp p 
h\2 
+( ) ™~p—sS SINa cos (@—pf 
p 
.\- sh (b\2 ” 
terms contamming p ( )>p > pl ) ete 17) 
p p p p 


ikeo iN a cos(@—B l(a do 


Pre p)U ik RY 


o\4 
terms containing ke ( ) > ke 
p 


({) ete. as 


Though this integral may be expressed by known 
functions for some current distribution functions 
gio), those funetions which will be used later 
in this report will, except in the case of \a 
constant, require rather complicated computations 
As ke | has been assumed, the term in the last 
bracket, containing Ae, may be neglected with a fair 
approximation. Thus, 


pret [ ©) — ite (£ ) y(a)do 


terms contalning ke (“) ete 19) 
p 


When using this approximation, the retardation 
along the antenna is neglected (but not from the 
antenna to the field point 
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hp 


Here, 
kh l, 


from which 
pPp-~ ‘ ihre : | ikea Sllli@cos (oO QO ( oO dla 
J, Cp) y 


terms containing } R ele 20) 


Using this approximation means to neglect the induc 
field term Also in this case it will gcenerally 
too cumbersome calculations to include thie 


tion 
involve 
second term in the square bracket; as in the cas 
with a fair approximation this 
This means neglecting th 


Then, 


considered above. 
term may be omitted 
retardation nloneg the antenna itself 


Poe — ere | (),) U(a)do 2 


In what follows, eq (16), (19), and (21) will be 
used as working formulas. The first and the third 
formulas ure both spec inl cuses of the second formula 
following work will be based) upon 
and specializations to the 
the end of this section 


treating 


Therefore, the 
formula (19 two other 
cases will be made only at 

\s was done by Wait (see 


a disk-loaded monopole, il Thies profitably be assumed 


footnote | In) 


here that the distribution funetion qd of the antenna 


current is expressed by the series 
2. A 


ly serting the above series: CXPression for the cur- 


rent distribution funetion  g(o in the integral 
expressing | a 
Pac DS A,p 


m= | [Cp) ie (ip) forte 


Vhis integral may be expressed by known functions 
value ol om The coethicients p 
a linear current distribu 


integral 
lor dealing with 
Lion, 1.e will be evaluated here. 

In order to carry out the integrations, .Y R¢ 


for owns 
HecessaLry 
p ana jp 
Is pul mto the standard form 

Y=C+2Re 
where 


COS a 


I Sillh aw COs (OO Oo 





Is Important to 


the integrations it 
/P?. It is found that 


In carrvine out 
know the sigh of N ih 


fp b\? /_p,.b 
V=(E) He) (ets) 
where 
yY=siIn a cos (o9—p 
COS @. 
As 


t mav easily be verified that 


Further defining 


the following is obtained 


lf hb i; ihe B+ i; 
) tan —tat 


© WL, VCJ wNL VN AY 


. ( tan e = tan a | 24 


\- \ 


These expressions are sufficient for calculating the 
contributions to the surface current density corre 
sponding to the terms contaming o” and o' in the 
current distribution g(e 


series expression for the 
ete corresponding to 


of the wire. The terms ps, ps. 
the higher order terms in the series expression lor 
gio) mav be found for example by using partial 
integration, but the calculation is tedious. In the 
following sections there is a use only for the terms p 
and p;, corresponding to a linear current distribution 
The calculation of Pis Ps etc., is therefore omitted 

Following the discussion given earlier, simplified 
expressions will be given here for po and p, applying 
when p is respectively very small and very large 
compared to the wavelength 


1 fi+B s \. 
di ; = + ah 
\ = \-\ \ f a 


hp 





lf b+ oe 
’ \ i , \ T \ ( ] 2 5b 
kp 
) hee ( s ) Poa 
. 4 = p ) Gb 
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3. Vertical Wire 


The main purpose of this paper is to calculate the 
earth current excited by a top loading consisting of 
one or more inclined or, as a special case, horizontal 
wires. However, a comparison of the contribution 
to the earth current from the top loading with the 
from the vertical member of the an- 
tenna The formulas for calculating this 
contribution may be obtained as a special case of 
the formulas developed above 

Considering un antenna of the height hand with a 
few top-loading wires of the length a, and assuming 
therefore a linear current distribution, the current J 
in the vertical member of the antenna may be ex- 
pressed by 


contribution 


is desired 


I—1,(1 —. 


current at the bottom of the vertical 
wire, and where xs In this seetion and the follow- 
ing sections J, will be reference current J 
and the heioht A of the antenna as the reference dis- 
tance d occurring in the above formulas. In using 
these formulas on the vertical member it is further 


where J, is the 


used as the 


set that 


The distribution funetion g(¢) for the current on the 


vertical wire is then expressed by 


) i 
ad he 
The following expression for the normalized earth 
current due to the vertical Wire is hereby obtained: 
Bite 
an h r. 
(). 
where 
Lp Li pye 
with 
A | 
| in 
| ; 27a 
ad 
| 4 
hi 
| heh | -_ 
y) tan 27b 


ih I+(7) 
yis(@) § > 


As the main purpose of this report is to investigate 
the ground current induced by the top loading, 
whereas the ground current induced by the vertical 
member is only calculated for comparison, it may be 
assumed, for the sake of simplicity, that the vertical 
wire carries a constant current 7; equal to the mean 
value of the actual current distribution, i.e 


Ar ed 
S09) 


IR 


We then have 


VW here 


and where Po IS expressed is above The final CX- 


pressions for the components of j then become 
F = | l | Nf | 

-7 p ' p\- p » 14 a 
Liv! () ) a ( i) 


] (). 2g 


tkh tan 


These expressions are equivalent to formula (7) in 
the paper by Wait (see footnote 1) referred to above 
In figure 6 |75) has been plotted for kA—0 as a 
function of ph and with a/h asa parameter. In this 
case, the last term in the first square bracket in the 
formula expressing } drops oul In the case of very 
large distances p, however, this last term may be- 
come significant and even dominant. Figure 7 
shows a similar plot of 7, for kk=—0.1. This value 
of kh corresponds approximately to the antenna for 
the planned U.S. Navy vif radio station, Cutler, 
Maine approximate data: f~15 ke, h~soo ft 


4. Top Loading 


The eround current excited yy an umbrella type 
of Lop loading Tih\ be calculated by adding vectort- 
ally the currents excited by the single wires In this 
section an expression is given for the surface current 
density J induced yy a single wire of the top lou ling 
Letting o=0 correspond to that end of the wire 
under consideration which ts connected to the verti- 
cal member of the antenna and o=1 to the outer end 
of the wire, the fact that the current at the end of the 
antenna is zero is expressed by 


g(1)—0. 
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) he ba ot th nl naa } ; norma ead noth aho } lop- oadin / é a 
ca 
1 paran 
Using the series expression introduced above for the In the case of a top loading consisting of one wire 
current in the antenna, it is found that if the wire is horizontal, an L-antenna), or of two 
wires (if the wires are horizontal, a T-antenna , the 
\ current on these wires will probably to a good 
> Ae =0 | 
— * approximation, be linear. Then, 
The current at o 0 is denoted by ae As in the last a 
section, the current J, at the bottom of the vertical l h 
wire will be used as the references current, J,. Thus : . " a 
| + 
h 
aio 


and all other coetherents 4 equal to zero 
On the other hand, if the top loading is of the 


In what follows, it will be asst “il that the 1 , 
ese stir imed or - oF load umbrella type, 1.e¢., 1f it consists of several equianhgu- 


Ing consists of m equiangularly spaced wires. Eacl 
‘er . 7 ’ larly spaced wires, the current on the wires will prob . 
wire will then earry the same current and conse bly } | > . P 
tly ably be better approximated by a quadratic function 
quen : 


In this cause, 


qi 


“~~ 
~ 


where 7, is the current at the top of the vertical wire h 
Using the series expansion for g(o 


and all other coeflicients A equal to zero The 


remaining part of this report will be concerned 
ii / mainly with the CUSseCS of il top loading consisting ol 
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1: Vv 
Jp 
6 
~ 
| IGURE GT Vorma eds we Cu ent densily Jo at the ground 
plane due to the ¢ ent nthe ertical mire mber of a fop- oaded 
antenna as a tunctior ot the normalized distance p h from 
the hase of the antenna and vith the normal ed Le nath al 
or the lop oadit / es as a paramete) the height h ot the 
antenna being given bukh 0.1 
' | . » ~ F | he curre 
Olle 7 or two (4 -) Wires and the current 
distribution will consequently be assumed to be 


linear. 
As Ink the last section, h Ww ill be used as the reference 
/ in ealeulating the normalized = surface 


distance d 
current densit Further setting 


the following is obtained from the formulas developed 


above 


| h h p > inn 
J a Silla cos (o—B)-+ cosa ll. 30a 
=~ fl ad ad 


BP. 


) SIN @ SID ré8) 30Ob 


where 
ad ad 
l h l h l I 
P } ; ‘ = 
7 a l f 7 a { \ A \ f 
h h 


p . h 
B =—— SIN a cos (QO B)—— cosa, 
a a 


For numerical calculations, it is of practical im- 
portance to have simplified expressions for P? valid 
when p is very small and very large compared to the 
wavelength respectively. Making use of the dis 
cussion given earlier, the following Is obtained 


kp l 
ad 
thd iB 
f 7 a NY’ ‘\ \ f v« i 
h 
] hp 
ad ad 
ae oe lh 7 
ion oa E Mie } an ( i) 
1+ 1-4 
h } 
ad 
ay] , 1 fk a 
ika (4) | etem Dy | a ka (~ ) ( > 
hi 


5. Horizontal Top Loading 


In practice a horizontal or nearly horizontal top 
loading is often used. In this section the simplifica- 
tion which mav be obtained in the formulas derived 
in the last section, as applied to the case of one or 
two wires when these wires are horizontal, will be 
considered, and a numerical calculation of the surface 
current will be carried out. 

Setting a=7/2 the following expressions are ob- 
tained for the components of the normalized surface 


current density 7“ induced by the horizontal wire 


¢ : é ) cos @—B).P sda 


2ria 


j= ("Yin @—8).P. 34h 


2r\a 


these formulas it is seeh, as has also been 


is parallel to the horizontal 


Krom 
shown before, that 7" 
wire. Therefore, j"=j" & may be 
unit vector pointing in the positive direction of the 
It is found that 


set, where s is a 


Wire. 
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The factor | Is expressed by the quantities RR. E. Vv. 
and VY, as in the last section. However. in the present 
cuse the eXpression for B reduces to 


The quantity mw 7 has been plotted In) figures Si 
and Sb for kA —0 and in figures 9a and 9b for kA =—0.1. 
and in both cases for ¢é—8—0° and 180° as a function 
of ph and with a/h as a parameter. From 
curves and from similar curves for other values of 
by multi 


these 


the p and the ® components of nj’ 


© De 
plying with cos (¢—8) and sin (@—8), respectively, 
are obtained 
r 
qa 
»/b 
pv 
ky Kt s Vi a / li 6? der } 
| pia / >t} ent lhe ho onta ( f 
; ; fy ” ho onta faop-loading u ‘ ! 1 ; 
fhe norma ed distance ph trom the base ot the rntlen? 
j i} this norma ead enol a? ol the fop-loadir 
‘ a ] 
5 B=1S 
6. L-Antenna 
An antenna of the heiht h and its horizontal 
member having the length @ as shown in figure 10 
are considered As in this case the number » of 


wires In the top loading is equal to 1, the contribu 


tion 7 to the normalized surface current densit\ 


coming from the horizontal top loading may be ob 
tained directly from figures 8 and 9 in the last section 
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1 
, 
- 
| I¢ REY Vo a / a ; lev a hie ound 
plane due tothe ¢ e; ; he ho ontal? pritee oranan ne 
eel one he onta lop-loading j ® as a tunctior nt the 
norma ed fistance pl iron fhe hase or the antenna and 
th the norma ed engih ah ot the lop-loading u easa 
paramete the heiaht h of the antenna ber Tiven f KI 0.1 
— a at 
’ 
vor re + LOK 
, zN Se KY Y/AN 
. Y 
PiGurRe 10 antenna 


The contribution ) obtained from these curves should 
be compared with the contribution J originating 
from the vertical member of the antenna as shown 
in figures 6 and 7. A comparison of these figures 
shows that, compared to the radial component 7°, of 
the normalized surface current density induced by 
the vertical the radial component 7" of the 
normalized surface current density induced by the 
horizontal vicinity of the 


Wire, 


wire is small also in the 








« antenna, and becomes insignificant at large distances. 
This was the same conclusion that was arrived at by 
Wait in his investigation of a disk- 
loaded monopole. A similar statement may be made 
regarding the azimuthal component 7/2 of the nor- 
malized surface current density induced by the hori- 
zontal wire. However, if the ground wires are 
placed radially around the antenna, as is usually 
the case, the eround losses induced by the azimuthal 
component 72 may be significant compared to the 
losses induced by even if . Is considerably smaller 
than 7 


(see footnote | 


In order to facilitate a comparison of 7: with 7 


the quotient 


a 
< sin (@—8 
j 

has been plotted in figures lla and 11b for kh=—0 

and for ¢—8—45° and 135 In the cases of d—B 

0° and 180°, respectively, the quotient is zero 

Curves have not been given for the case of kh=—0.1 

treated earlier, as these curves would deviate only 

slightly from the curves in the case of kh=0 over 

the significant range of values of p/h. 

7. |-Antenna 
\ T-antenna of the height h and with its hori- 
zontal member having the leneth 2a as shown in 
és figure 12 is considered. Letting 8=8, denote the 
angular position of one of the two horizontal wires 
ot which the top loading consists, we have for the 
— 
~ — ~ >, a 
4 
So 
“ 
¢ 
a 
Q= 3. 
* 
. 
~ 
" 

FIGURE 11 Re veen the azimuthal component of 
the ace ; sity due to the current in the horizontal 
mie mile ) ! antenna and fhe face ‘ ent densil 

fide 0 ; 1 membe ‘ / same antenna, 
It hi I | 
he ! B=45 4 
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the other radial wire 


angular position 6= 6, of 
Utilizing the symmetry properties of the 


b2= 6 T. 
normalized surface eurrent density )' indueed by 
one radial wire it is then found that the following 
expressions for the normalized surface current 
density 7” of the total top loading by adding the con- 
tributions from each of the two wires, 


where s,; and Ss, denote unit vectors in the positive 
directions of the two top-loading wires, respectively, 
as shown in figure 12, and where 7? (w) is the funetion 
introduced in the section on a general, horizontal top 
loading. Utilizing 





and 
} ij ds WV 
it is found that 
14 (gd 3 S , o 3 va 
where 
Is, (@—-8 Js (6-8 Js (7>—(@—BP s 
= a ann aan Q ae 
<—_— — 
| A A 
Se S| 
h 








WKLNY IDS ASP 


FiGguRi 


LS Ler hh whiwlis 


— 
12. |-anlenna 


As the number of wires in the top loading is equal to 
a the funetion 


mav be obtained as one-half the value given by the 
curves in figure 8 or 9. 


It is easily verified that dl 8) satisfies the fol- 


lowing symmetry relations, 


9 
(54a 


9 
(39b 








It follows from the above expression for 7” that 


¥ T 
j’( +5 )=0 


Because of the svmmetry relations mentioned abov e, 


it is necessary only to calculate the normalized 
surface current density for @— 8, being in the interval 
Oto wr? 


The quantity |) 7"%(@—8) | has been plotted in 
13 for kh 0 and for o D 0° as a function 
and with a/h as a parameter. This diagram 
giving the contributions to the normalized surface 
current density arising from the horizontal wire 
should be compared to the contributions from the 
vertical member shown in figures 6 and 7. This 
comparison leads to the same qualitative statements 
regarding the significance of the relative contribu- 
tions from the horizontal wire as were arrived at in 
L-antenna 


figure 
of ph 


the case of the 


hiGguRt 13 Vormalized surtace ¢ ent density , at the 
yround plane due to the current in the ho ontal membe ot 

a antenna as a tunction ot the norma ed distance p} 
om the base ot the antenna and with the norma ed lengli 


ah of the top-loading wires as a paramete 


It med th +) } hr bof th nienna nfinit 
‘ ‘ he enna ely ( | 


In order to be better able to survey the situation 
in figure 14 the following quotient is plotted, 


for kh —0O and for ¢6=—45° as a function of p A and with 
ajh asa parameter. For @—0° and 90° the quotient 
Is Zero Pay Ing attention to this facet and comparing 
figure 14 with figures lla and IIb. it can probably 
be said that in general the azimuthal component of 
the surface current density induced by the top 
loading is smaller for a antenna with a horizontal 
wire of the length 2a than foran L-antenna of the same 
heiht and with 2 horizontal wire of the leneth a 


tha 


| IGURE 14 Ratio between the azim component or the 


surface current 


de nsil 


members of a leaqnienna and the surface cH ent 


j due to the ertical membe ol the same ante 


On the other hand, the radial component will at 
some points be larger, at some points smaller for the 
T-antenna than for the corresponding L-antenna 
For an antenna with a radial ground wire system i 
is essential to keep the azimuthal component of the 
surface current density as small as possible in orde1 
to make the total cround losses the smallest possible 
In this respect a T-antenna is probably in general 
preferable to an L-antenna 


8. Antenna with Four Horizontal 
Top-Loading Wires 


The theory developed above for the surface cur- 
rent induced on the ground plane by a straight wire 
with a linear current distribution will be used also 
for finding the earth currents induced by a top- 
loaded antenna with four equiangularly spaced 
wires in the top loading. Although in the present 
case the assumption of a linear current distribution 
aus Was mentioned before, probably does not aupprox- 
distribution so well as in 
Is desirable Lo 


imate the actual current 

the case of an T-antenna, it 

include the case of an antenna with four horizontal 

top-loading wires, as this case is very close to the 

case of a disk-loaded antenna; on the basis of the 

should therefore 

between the wd 

the theory 
However. 


L- OF a 


this section it 
comparison 
report ana 


results obtained in 
be possible Lo make 2 
theory developed in this 
developed for a disk-loaded monopole 
it should be kept in mind that whereas the formulas 
and numerical results obtained in this report are 
based on the assumption of a linear current distribu- 
tion on the top loading, Wait in his investigation of 
a disk-loaded antenha assumes a quadratic eurrent 
distribution on the disk 

The length of the vertical member of the antenna 
is denoted by kh and the leneth of the top-loading 
wires by a as indicated in figure 15. Further, the 
azimuth of the top-loading wires is denoted by 8 


B., Bs, and By, respectively, and the unit’ veetors 
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It may be verified that 7*(@—,) satisfies the follow- 
ing svmmetry relations: 





I (—(@— Bi S1Js (@— Bi) — 8295 (@— 2B), ta 
V( ; (@—p ) . J o— 2, S. 1, (@—B,;), 12b 
7 r™—(9—O Sis (@— 8 Soh, (O—B {2¢ 

It follows from the above expression for) that 


(F )} ib 


where 6 denotes a unit vector pointing in the radial 
n-loadina wires direction. As the number » of wires in the top load- 
ing is equal to 4, the function 








pomting im the positive direction of these wires by )}.(¢—6 
S., Sg, and s Then, , 
may be obtained as one-fourth the value given by 
- the curves in figure 8 or 9. 
nt =) Because of the above relations the whole current 


field will be known, if it has been caleulated in the 
angular space 


) y) T - 
() o 5 * 
' 1 
) e) _ 
) , ° . 
. Using the svmmetry relations stated above it is 
ind found that » (@—3,) is purely radial for 
‘ ‘ . 
oO ) )—» 
4 
where Pp is an integer For these angles the radial 
Making use of the sVinmetry properties of the nor- component 
malized surface current density 7” induced by one i+(@—B 


radial wire it is found that the following expressions 
, = is given by 
for the normalized surface current density J of the ; 


total top loading 





2G )-a(F ) J for 6 B a t45 


where q is an integer. $4 


A 
where 


The quantity 7, (¢@—8,) has been plotted in figure 

16 for kh~O and for @—8,—0° as a function of p/h 

o j ( - o—8 ) ; ( ¢ . ) tb and with a/h as a parameter. These curves should 
- 2 be compared with the curves given in Wait’s report 
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(see footnote 1) on the earth currents induced by the 
disk in a disk-loaded monopole. It is seen that the 
currents in figure 16 have the same trend and the 
same order of magnitude as the curves computed by 


Wait. 


p 
hy Kk lb Radial con ponent of norma ed surtace 
ni density due to the ¢ ent ; fhe ho onta mer 
Qnlenna will oO ho onta fop-loading (SS as a lune 
ol the norma ed distance pt or hel eolthe ar 
} } a j i} ‘ / oO ; mp lil 





9. Conclusion 


Based on a formula for and numerical calculations 
of the surface current density on a ground plane 
above which is placed a straight piece of wire with 
a length that is small compared to the wavelength 
and with a linear current distribution, the contribu- 
tion to the ground current of the top loading of 
L- and T-antennas that are small compared to the 
wavelength has been found and plotted Whereas 
the contribution from the top-loading to the radial 
component of the surface current density, as was 
shown earlier by Wait, in general may be neglected, 
as it is insignificant compared to the contribution 
from the vertical wire, the azimuthal component of 
the surface current density produced by the top 
loading may be important in calculating the ground 
losses in the case of a svstem of radial ground wires, 
as it is perpendicular to these wires. ‘The formulas 
and numerical results in this report form a basis for 
calculating the eround losses due to this azimuthal 
component 


This investigation was carried out while the author 
was temporarily emploved as a consultant at the 
National Bureau of Standards in’ Boulder. He 
thanks the Bureau for the Opporvunity to do this 
work He also thanks James R Wait for sugevesting 
the problem and \Iaxine W Clark fol 


numerical calculations 


making the 


BoutpeEr. Co.o.. Januarv 27. 1959 
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Important Notice 


Beginning in July 1959. the Journal of Research will be published in four 
separate sections which may be subscribed for individually. This change is being 
made to provide a publication whose subject matter concentrates on particular 
fields of interest, thus meeting the specialized needs of individual scientists. engi- 
neers. and mathematicians. 

Section A. Physics and Chemistry, to be issued bimonthly. will present 
papers of interest primarily to scientists working in these fields. This section will 
cover a broad range of physical and chemical research, with major emphasis on 
standards of physical measurement, fundamental constants. and properties of 
matter. 

Section B. Mathematics and Mathematical Physics, to be issued quar- 
terly. will present results in pure and applied mathematics. including mathematical 
statistics, theory of experiment design. numerical analysis, and short numerical 
tables: theoretical physics, chemistry. and engineering. where the emphasis is on 
the mathematical content or methodology: and logical design, programing. and 
applications of electronic computers and computer systems. 

Section C. Engineering and Instrumentation, to be issued quarterly. 
will report results of interest chiefly to the engineer or applied scientist. This sec- 
tion will include many of the new developments in instrumentation resulting from 
the Bureau’s work in physical measurement, data processing. and the development 
of test methods. It will also cover some of the work in acoustics. applied me- 
chanics, building research, and cryogenic engineering. 

Section D. Radio Propagation, to be issued bimonthly. will report research 
in radio propagation, communications. and upper atmospheric physics. ‘Topics to 
be covered include propagation in ionized media. scattering by turbulence, the 
effect of irregular terrain on propagation. diffraction and scattering by solid 


obstacles. propagation through time-varying media. and antennas. 


All present subscribers to the Journal will receive an announcement from the 
Superintendent of Documents when advance orders for the new sections can he 


accepted, 
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